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INVESTIGATION OF THE FORCES THAT GOVERN THE THREE-DIMENSIONAL 
PROPAGATION AND EXPANSION OF CORONAL MASS EJECTIONS FROM SUN 
TO EARTH 

Robin C. Colaninno, PhD 
George Mason University, 2012 
Dissertation Director: Dr. Angelos Vourlidas 

In the last few decades, we have discovered that the environment of our solar system is 
as dynamic as terrestrial weather. The source of space weather is the Sun which produces 
winds and storms that effect modern human systems. The most geo-effective aspect of 
space weather is Coronal Mass Ejections (CMEs) which are analogous to terrestrial hurri- 
canes. These powerful storms, comprised of plasma and magnetic fields, can significantly 
disrupt Earth's magnetic field and cause a range of terrestrial effects from the aurora to 
the destruction of technological infrastructure. It is only recently that we have been able to 
continuously monitor CMEs as they progress from the Sun to Earth. Even with continuous 
monitoring with remote sensing observations, we are still unable to accurately predict the 
arrival or terrestrial impact of a CME. There is much about the evolution of CMEs we do 
not understand. 

In this study, we analyze nine CMEs from the Sun to Earth as observed in both the 
remote sensing and in situ data sets. To date, this is the largest study of Earth impacting 
CMEs using the multi-view point remote sensing and in situ data. However, the remote 
sensing and in situ data of the same CME cannot be directly compared. Thus, we use 



several models to parameterize the two data sets. We are able to compare the arrival 
time, Earth impact speed, internal magnetic field, size and orientation as derived from the 
remote sensing and in situ methods. By comparing these results, we hope to gain a more 
comprehensive understanding of the inner heliospheric evolution of CMEs. 

In this study, we track CMEs from the Sun to 70% - 98% of the distance to Earth 
with the remote sensing data. We analyze the propagation and expansion of the CMEs by 
applying two geometric models of their structure. From the derived kinematics, we compare 
the predicted arrival times and impact velocities with the in situ data. We find that even 
with nearly continuous observations and the best available model of the CME structure, 
there is still a significant error in the predicted values. We discuss the possible causes of 
these errors. 

To investigate the drivers of the CME, we use the derived kinematics and mass of the 
CMEs to calculate the net force. We estimate the various forces acting on the CME as 
predicted by three theoretical models of CME propagation and expansion and compare 



these results with the observational results. We find that the flux rope model of Chen 



(1989) provides the best agreement with the observations. With the flux rope model, we 
are able to predict the internal magnetic field of the CME near Earth from the remote 
sensing data to an order of magnitude. This result is of great importance to space weather 
predictions since the strength of the internal magnetic field is a major factor in the geo- 
effectiveness of a CME. We discuss a possible method of better fitting the flux rope model 
to remote sensing observations with the goal of improving the accuracy of the predicted 
magnetic field values. 

Finally, we compare the size and orientation of the CMEs as predicted from the remote 
sensing and in situ data. We find very little agreement between the values derived from the 
two data sets. This aspect of the study also has large implications to space weather since 
the orientation of the CME's magnetic field is the other major factor in geo-effectiveness. 
Since the in situ data is our only source of direct information of the magnetic field and 
remote sensing data the main means of prediction, it is essential that the analysis of these 
two data sets is reconciled. 



Chapter 1: Introduction 



Coronal mass ejections (CMEs) are complex, dynamic, three-dimensional structures com- 
prised of energetic particles, plasma and magnetic field that are ejected away from the Sun. 
In coronagraphic and heliospheric image data, CMEs appear as a structured brightening in 
the corona moving away from the Sun. Many large CMEs have a three-part structure which 
includes a bright leading edge, a dark void, and a bright complex core, see figure [lT] ( [Illing 



and Hundhausen 1985). The bright relatively smooth leading edge is created by coronal 



material displaced by the motion of the CME. The leading edge demarcates the boundary 
of the magnetically closed region of the CME. A shock is sometimes associated with the 
front of the CME. The shock forms ahead of the CME leading edge and is much dimmer 
than the leading edge. Near the center of the CME envelope is a bright core often with 
a complex, twisted structure. The CME core is associated with material from an eruptive 
prominence. Between the leading edge and the core is a dark region comprised of rarefied 
material suggesting a magnetically closed cavity. 

The kinematics of CMEs is driven by the forces exerted by the magnetic field and 
gas pressure both internal and external to the CME. Internal to the CME, its magnetic 
structure and gas pressure drive the motion of the CME. These internal forces interact with 
each other and with the external magnetic field and gas of the solar wind. As a CME leaves 
the Sun, its evolution can be separated into two components; propagation and expansion. 
The propagation speed is the change in height of the CME's center of motion above the 
solar limb. The expansion motion is the movement of the leading edge way from the center 



of motion (Colaninno and Vourlidas 2006). Figure 1.2 shows a diagram of the change in 
height and radius of the CME as it moves in the coronagraph's field of view. A parameter 
often used to quantify CME motion is the speed of the leading edge. Ideally, the leading 
edge speed is the maximum of the superposition of the propagation and expansion speeds. 
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Figure 1.1: A three-part structure CME showing the loop hke leading edge, dark cavity and 
bright core as seen by SOHO-LASCO C2 on December 2, 2002. 



Analytical models can be used to describe how the magnetic field and gas pressure drive 



the propagation and expansion velocity of CMEs (Chen 1989, Low 1982a, Pneuman 1984 



Siscoe et al. 2006). 



Near Earth, magnetic and gas parameters can be detected by in situ instruments. Several 
spacecraft have instruments that detect the electromagnetic field, particle composition, 
speed, pressure, temperature and density of the near Earth space environment. Signatures 
in the interplanetary medium have been linked to the passage of the spacecraft through 
a CME. These signatures seen in the in situ data are referred to as interplanetary CMEs 
(ICMEs). Some ICMEs have a particular well-organized structure called a magnetic cloud 
(MC). The signature of a MC is defined by enhanced magnetic fields, a smooth rotation of 
the magnetic field vector, and a low proton density and temperature ( Burlaga et al.||T981 ). 

The magnitude and direction of the magnetic field of the ICME as it impinges on Earth's 
magnetosphere is the main driver of space weather and geomagnetic storms. Several studies 



found that large geomagnetic storms are caused by ICMEs (Gosling et al. 1991, Richardson 



et al. 2002). However, due to the lack of observational data and the complex nature of 



2 









?) 












\ w / 




R / 


^ t ^ 


^ t+dt ^ 


t+2dt ^ 



Figure 1.2: As the CME moves away from the Sun, the center of motion (CM) increases 
in height (R). The change in R with time is the propagation speed. At the same time the 
distance between the CM and leading edge increases (r). The change in the CME radius, 
r, is the expansion speed. 



interplanetary interactions, the development of CMEs into ICMEs is unclear. Based on 
observations between 1996 and 2000, unique CME-ICME associations have been made for 



approximately half off all ICME events (Zhang et al. 2003). 



1.1 Observations of CMEs 

In this study, we use two types of data to study CMEs in the corona and heliosphere; 
remote sensing and in situ data. Remote sensing observations are made by telescopes 
that image solar light scattered by electrons in the corona. In 1972, the first CME was 



observed in remote sensing data (Tousey and Koomen 1972) with the coronagraph aboard 



the Orbiting Solar Observatory-7 (Koomen et al. 1974). Since then, CMEs have been 



observed by a number of space-borne coronagraphs. Here, we use remote sensing data from 
two spacecraft; Solar TErrestrial RElations Observatory (STEREO; Kaiser et al.||2008 ) and 



SOlar and Heliospheric Observatory (SOHO; Domingo et al. 1995). In situ observations 



are made by instruments in the heliosphere measuring properties of the corona plasma as it 



passes over them. Burlaga et al. ( 1982 ) was the first to associate a magnetic cloud measured 



in situ by Helios 1 spacecraft with a CME seen in remote sensing data with the Solwind 



coronagraph on the spacecraft P78-1 Sheeley et al. (1980). Here, we use in situ data from 



Figure 1.3: Shows LASCO C2 (left) and LASCO C3 (right) observations of a large CME on 
16 August 2001. The LASCO C3 image also shows the impact of fast, high energy protons 
that struck the SOHO spacecraft's detectors. 



the WIND mission to find the corona plasma signatures of CMEs near Earth. 
1.1.1 Remote Sensing 

Since its launch in 1996, the Large Angle Spectroscopic Coronagraph (LASCO) aboard the 
SOHO mission has been observing the corona, almost continually. LASCO observes the 
corona from the L-1 Lagrange point along the Sun-Earth line. LASCO is comprised of two 
operating coronagraphs; C2 and C3. The LASCO C2 has a field of view from 2.2 to 6 Rq. 



LASCO C3 has an overlapping field of view from 3.8 to 32 Rq. Figure 1.3 shows an example 
of the LASCO C2 and C3 observations during a large CME. 

The Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI) instru- 
ment suite abroad the STEREO mission offers a novel view of the Sun-Earth system. 
STEREO, launched on October 25, 2006, is comprised of two spacecraft with nearly identi- 
cal instrumentation. The two STEREO spacecraft orbit the Sun; the Ahead (A) spacecraft 
slightly faster than Earth and the Behind (B) spacecraft slightly slower. The two spacecraft 
are separating from Earth at a rate of 22.5° per year. SECCHI is the imaging instru- 
ment suite of the mission. The SECCHI suite includes an extreme-ultraviolet imager, two 
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Figure 1.4: Shows SECCHI instrument suite's field of view along the ecliptic plane from 
the Sun to Earth for a single point of view. The insert shows the EUVI, CORl and C0R2 
fields of view which cover the corona from the solar disk to 15 Rq approximately equal to 
4° elongation. The HI-1 field of view begins at 4° and is 20° by 20°. The HI-2 field of view 
is 90° square and a slightly overlaps the HI-1 field of view. The field of view of both HI 
instruments are centered on the ecliptic plane. The trapezoidal occulter at the end of the 
HI-2 field of view is to block Earth shine. 



coronagraphs and two heliospheric imagers. Together these instruments image the inner 
heliosphere from the Sun to Earth from two viewpoints. The two coronagraphs in the SEC- 
CHI suite are CORl with a field of view from 1.5 to 4 i?Q and C0R2 from 2.5 to 15 Rq 



(Howard et al. 2008). SECCHI also includes two heliospheric imagers (HI-1, HI-2) which 
are similar to coronagraphs but have no occulter and a field of view off pointed from the 
center of the Sun. The heliospheric imagers view one side of the corona along the Sun-Earth 
line. HI-1 and HI-2 have square fields of view centered on the elliptic plane from 15 to 84 
Rq (20°) and 66 Rq to 1 AU (70°), respectively ( [Howard et al][2008l ). Figure [O] shows the 
combined fields of view of the SECCHI instrument suite. 

When studying CMEs with remote sensing data it is very advantageous to combine 
the data from SOHO and STEREO. The LASCO data set has offered abundant insight 
to the nature of CMEs over the solar cycle. From the LASCO data, physical parameters 
such as position, angular width, speed, mass and acceleration have been cataloged in the 



CDAW database for all observed CMEs ( [Yashiro et al.||2004D . However, the LASCO CME 
measurements were made using data from a single view point and are two-dimensional 
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projections of three-dimensional values. With the two views of SECCHI, we are able to 



derive the true three-dimensional position, velocity and mass of CMEs (Colaninno and 



Vourlidas 


2009 


Mierla et al. 


2010 



have not been measured continually from the Sun to 1 AU. Thus the propagation of CMEs 
in the inner heliosphere and three-dimensional structure and dynamics of CMEs has only 
become available with the SECCHI data. For the first time, we can observed CMEs from 
the Sun to 1 AU from two viewpoints away from the Sun-Earth line. 



1.1.2 In Situ 

The WIND spacecraft was launched on 1 November 1994 with the primary science objectives 
to provide data for magnetospheric and ionospheric studies and investigate basic plasma 
processes occurring in the near-Earth solar wind. The WIND spacecraft has a suite of 
several instruments to collect data in situ. In this study, we will use data from two of 
these instruments; Magnetic Field Investigation (MFI; Lepping et al.|[l995 ) and Solar Wind 



Experiment (SWE; Ogilvie et al. 1995). The MFI instrument is a triaxial magnetometer 
which provides the magnitude and direction of the solar wind's magnetic field. The MFI 
instrument has a dynamic range of ±4 nT to it 65,536 nT with a sensitivity of 0.008 nT/step 
quantization. The SWE instrument includes two Faraday cup sensors and an electrostatic 
analyzer. The Faraday cup sensors provide the density, velocity and temperature of the ions 
of the solar wind. The electrostatic analyzer measures the solar wind electron pitch-angle 
distribution function. With the data from these two instruments, we can detect CMEs and 
MCs as they pass over the WIND spacecraft. 



1.2 Theoretical Models of CME Structures 

In this study, we use three theoretical models to analyze the three-dimensional structure of 
CMEs from the observational data. We use two models to fit the remote sensing data; the 
Elliptical Model and the Graduated Cylindrical Shell Model. To model the in situ data, 
we use the Force-Free Elliptical Flux Rope Model. We have chosen to use the Elliptical 
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Model because the remote sensing observations show the cross-section of the CME in the 
inner hehosphere to be elhptical. For the same reason, we have chosen an in situ model 
which also has a elliptic cross section. We expect that as the CME expands latitudinally 
the trajectory of an in situ spacecraft through the CME could be approximated by an 
elliptical cylinder. However, the Elliptical Model is very limited especially in dealing with 
the projection effects in the remote sensing data. Thus, we also use the more sophisticated 
Graduated Cylindrical Shell Model to determine the trajectory of the CME. The drawback 
of the Graduated Cylindrical Shell Model is it has a circular cross-section. Thus, we can 
not use this model to analyze the expansion of the CME in different directions. 

1.2.1 Elliptical Model 

The most common method of deriving physical parameters from remote sensing images is 
direct measurements made of features in the image. These measurements are made by an 
observer choosing a feature in the image and tracking that feature over many images. The 



data from the SOHO LASCO CME Catalog ( [Yashiro et al]|2004| ) is collected in this way 
where the speed, acceleration, position angle, and width are determined by an observer 
selecting the pixel locations of the leading edge height and angular extent of the CME over 
three or more images. 

The structure of the CME can also be analyzed using direct measurements by selecting 
the points in the image which define the envelope of the CME and fitting a geometric 
structure to these points. The geometric structure most widely used in the literature is 



the Cone Model. Howard et al. ( 1982 ) introduced the Cone Model and it has been widely 



Xie et al. 


2006 


Zhao et al. 


2002 



Model subscribes an ice cream cone like structure to CMEs where the cavity of the CME is 
described by a sphere and is attached back to the surface of the Sun by a spherical cone. The 
geometry of the Cone Model is represented in figure [L5) [a). The Ellipse Model used here 
is a two dimensional variation of the Ellipsoid Model were the cavity of CME is described 
with an elliptical cross-section. Figure |1.5[b) shows the geometry of the Ellipse Model. Of 



a) b) 

Figure 1.5: Geometric representation of the (a) Cone Model and the (b) Elhpse Model. For 
the Cone Model the radii of all the axes are equal. For the Ellipse Model the three axes 
can be of different lengths. 



course this model is easily reduced to the Cone Model where the length of all the axes are 
equal. 

The Elliptical Model is implemented by making direct measurements of the outer en- 
velope of the CME structure in the image. The trailing (sunward) edge of the CME is 
the most difficult to identify. We define the trailing edge as the apex of the V feature of 



the CME (Kahler and Webb 2007). If there is no clear V feature, we identify any trailing 
feature that is easily track in the images. After six points have been selected, an ellipse of 
the form equation |1.1| is fitted to the direct measurements and over plotted on the image. 
The observer can then acept the fit or add more points until a satisfactory fit is obtained. 
Thus the number of points used for each ellipse fit is determined by the observer. 

x{6) = Xc + a cos 6 cos cj) — b sin 9 sin cj) 

y{9) = yc + acosOsmcf) + hsinO coscj) (1-1) 



For each fit, we get the center of the structure, {xc,yc), the radii of the major and minor 
axis, (a, 6) and the tilt angle {(f)). We define the minor axis as the parallel radius and the 
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major axis and the perpendicular radius where parallel and perpendicular are defined with 
respect to direction of the CME propagation. 

The Elliptical model is a very basic first order approximation of the three-dimensional 
geometry of a CME. While the simplicity of the model is compelling, it is merely a geometric 
construct of the density and has no theoretical explanation. Also the model does not address 
the effects of projection in the images since it relies on direct measures. Direct measurements 
of CMEs are useful because of their simplicity of collection. However, direct measurements 
must be interpreted with care because of the projection effects that are dependent on the 
orientation of the CME. All measurements made directly from remote sensing data are 
projected on the plane of the observation. For example, the perpendicular diameter from 
the direct measurements is an upper limit to the cross-section of the CME. If the CME 
is oriented face-on then the perpendicular diameter is the width of the CME. Only if the 
CME is oriented edge-on is the perpendicular diameter the cross-section of the CME. For 
any orientation of the CME between these two extremes, the perpendicular diameter will 
be larger than the cross-section. It is difficult to correct for this enlargement since, we do 
not know which part of the CME will be projected onto the plane of the sky as the furthest 
extent of the flank. A similar argument can be made about the parallel diameter from 
the direct measurements. However, in contrast to the perpendicular diameter, the parallel 
diameter is the lower bound of the CME cross-section. If the CME is oriented in the plane 
of the image then the parallel diameter will be the cross-section. If the CME is oriented 
out of the plane of the image then the parallel diameter will be shortened by the projection 
onto the plane of the image. Thus the parallel diameter will always be less then or equal 
to the true cross-section. 

1.2.2 Graduated Cylindrical Shell Model 



The graduated cylindrical shell model (CCS) was developed by Thernisien et al. (2006) and 



Thernisien et al. (2009). This model provides a method for analysis of the three-dimensional 



morphology, position and kinematics of CMEs in white-light remote sensing observations. 
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CME Model 
Skeleton - 





Figure 1.6: Representation of the graduated cylindrical shell (GCS) model (a) face-on and 
(b) edge on. The dash-dotted line is the axis through the center of the shell, C. The solid 
line represents a planar cut through the cylindrical shell at the origin. The width, W, of 
the model is defined as the largest vertical extent, dotted line. The radius, a, defines the 
circular cross-section. These values are controlled by the model fitting parameters height 
(H), ratio (5) and half angle (a) ( Thernisien et al.|2009 ). 



The GCS model uses forward-modeling techniques which allows the user to fit a geometric 
repesentation of a flux rope to CME observations. The geometry of the empirically flux rope 



model is represented in flgure 1.6 The model is positioned using the longitude, latitude 
and the rotation parameters. The origin of the model is fixed at the center of the Sun. 
The size of the flux rope model is controlled using three parameters which define the apex 
height, foot point separation and the radius of the shell. The details of the model as well as 
the derivation of many secondary parameters used in this study are discussed by IThernisien 



(2011). 



The GCS model has been used in several studies of various aspects of the three- 



dimensional structure and kinematics of CME in the heliosphere. Lynch et al. (2010) 



successfully used the GCS model to derive the orientation of the flux rope in remote sensing 



data for comparison with in situ observation. Vourlidas et al. (2011) used this model to 



show the rapid rotation of a CME in the middle corona. Also, Poomvises (2010) used the 
GCS model to analyze the propagation and expansion of several CMEs with the goal of 
deriving a solar wind drag force. 
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This model is implemented by over plotting the projection of the cylindrical shell struc- 
ture onto each image. The observer then adjusts the six parameters of the model until a 
best fit with the data is achieved. We fit the GCS model to all available images from the 
CME's emergence in the SECCHI C0R2 and LASCO C2 fields of view until the SECCHI 
HI-2 images. We analyze the variation of the model fit parameters over time to derive the 
kinematics of the CME. 

The GCS model is a more sophisticated method of modeling the three-dimensional 
structure of the CME when compared with the Ellipse model. The GCS model correctly 
handles the effects of projection for each image's point of view. Also the geometry of the 
GCS model is a good proxy for the fiux rope like magnetic structure of CMEs. However, this 
technique is limited in that the empirically defined model has a circular cross-section. Thus 
we are only able to fit the overall structure of the CME and cannot model any distortion 
in the interplanetary space. 

1.2.3 Non Force-Free Elliptical Flux Rope Model 



Since Burlaga et al. (1981) identified the signatures of magnetic clouds, several models 



have been presented in the literature for reconstructing their magnetic structure. The most 



common is the model of Lepping et al. ( 1990 ) which solves a force free-form of Maxwell's 



equations in cylindrical coordinates. The force-free cylindrical torus is the basic framework 
for understanding the interplanetary magnetic cloud. It has provided an approach for the 
interpretation for the magnetic topology and orientation of the magnetic cloud from in situ 
data. 



The analytical model used in this study with the in situ data was published by Hidalgo 



et al. (2002) and developed by Nieves-Chinchilla et al. (2009) with successive improvements 



in order to facilitate implementation. It is similar to the Lepping et al. (1990) with two 



changes that make it more general: 1) a non force- free condition is assumed; and, 2) the 
coordinate system chosen to resolve the Maxwell equations is elliptical cylindrical. Both of 
these conditions have consequences in the physical understanding of the magnetic cloud. 
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The first condition has consequences in the overall picture of the CME evolution in the 
interplanetary medium. Locally, at 1 AU, the assumption that the system is under a force- 
free condition is not correct. The second condition has a direct impact on the geometry of 
the magnetic cloud. The elliptic cylindrical coordinate system allows the magnetic cloud 
to have an elliptic cross-section as indicated by the remote sensing data. The relation of 
elliptical cylindrical coordinates to cartesian is given by 

X = r cosh 77 cos (p 

y = y 

z = r sinh ?7 sin (1-2) 

where r is the radial coordinate, ip is the angular coordinate where 99 G [0, 27r) and ij is the 
coordinate that controls the eccentricity of the cross-section. 

The non force- free Elliptical Flux Rope model can generate many parameters to describe 
the CME. In this study, we are focused on comparing the orientation and size of the magnetic 
cloud with the remote sensing data. Prom the fitting of the model to the experimental data, 
we get several parameters related to the orientation of the flux rope. The three angles, {9, 
(f), ^) are the Euler angles relating the Y, X and Z axis, respectively, of the magnetic 
cloud to Geocentric Solar Ecliptic (GSE) coordinates. Also the vector, gives the impact 
parameter, distance of the spacecraft crossing from the center of the flux rope. The shape 
of the flux rope is given by two parameters, radius of the major axis and the ratio of the 
major and minor axis given as a percentage. 

1.3 Theoretical Models of CME Propagation and Expansion 

Even with the availability of data from remote sensing and in situ instruments, without the 
application of analytical models, we still do not have access to the underlying physics of 
CME dynamics. Many theoretical models of CMEs allow us to relate their propagation and 
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(a) 



(b) 



Figure 1.7: The relation between the coordinates system of the magnetic cloud 



^NM ^NM r^NM^ and the GSE coordinates system {Jp^^ ^^^^ ,Z^^^). The angles given 



by the model are (a) the longitude (^), (b) the latitude (^) and (c) the tilt angle (^) 



expansion back to unobservable quantities of the CME such as the mass, dimensions and 
kinematics. The three propagation and expansion models we will use here are; kinematic, 
self-similar MHD and flux rope. Each of these models takes a very different approach to 
solving for the propagation and expansion of CMEs. The models also have varying levels 
of complexity and physical insight. 

1.3.1 Basic Equations of Magnetohydrodynamics 

To begin reviewing the theoretical models of propagation and expansion; we must first 
review the basic equations of magnetohydrodynamics. The movement of a CME is controlled 
by the balance of forces interior and exterior to the CME. The two main forces that are 
interacting in the solar corona are the gas pressure and magnetic field forces. Due to the 
high coronal temperatures, the corona is completely ionized. The corona is a magnetized 
plasma and can be described by magnetohydrodynamics (MHD) to the first order. MHD 
is the combination of Maxwell's equations. Ohm's law, and the conservation equations of 
fluid dynamics to describe the bulk properties of a one-fluid plasma. Ideal MHD assumes 
that the plasma has infinite conductivity (or no resistivity). This assumption allows us to 
eliminate the electric field and current density terms from Maxwell's equations. Thus the 
plasma is parametrized by its density, pressure, magnetic field and velocity. 

All of the analytical models reviewed here, to some extent draw on the basic equations 



( [Nieves-Chinchilla et al.|[2005| . 
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of ideal MHD. Here we will briefly review the basic equations of MHD for later use in the 
CME propagation and expansion models. The self-similar MHD model solves the ideal MHD 
equations directly while the kinematic and flux rope models use the same simplification of 
Maxwell's equations used to derive the ideal MHD equations. 

Maxwell's equations describe the evolution of the electric field, E(r, t), and magnetic 
field, B(r, t) in response to the current density J(r, t) and the electric change density g(r, t): 



Faraday's induction equation : 

Ampere's law : 
Gauss's law: 



V X E 



1 dB 

c~dt 



V X B = 47r J + 



V • E = 47rg 



c~dt 



(1.3) 

(1.4) 
(1.5) 



Gauss's law for magnetism : 



V • B = 0. 



(1.6) 



To relate the motion of the plasma with Maxwell's equations, we also need to include Ohm's 
law 

E+ ^(v X B) = ??J. (1.7) 

c 

where rj is the electrical resistivity and v is the mean velocity of the plasma. The equations 
of gas dynamics for the the evolution of the density, p(r, t), and the pressure, p{r, t) are: 

dp 

continuity equation : ^ ^ ^^^^ ~^ ^^'^"^ 

/(9v \ 

equation of motion : i W ^ / ~ + J x B + pg (1.9) 

equation of energy : dt \~ ) ~ ^ (-'-•-'- ^) 



where g is the acceleration due to gravity and 7 is the adiabatic index. 

However, these equations are not yet in the most convenient and compact form for the 
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self-consistent study of plasma. We begin reducing the equations by assuming the ideal 
case where the resistivity in equation L7 is r/ = 0. Thus, we can substitute equation L7 
into equation 1.3 to eliminate the electric field term, 



dB 

~dt 



V X (v X B). 



(1.11) 



We can easily assume that the plasma is non-relativistic, f <C c. We can then estimate the 



order of magnitude the electric and magnetic field terms in equation 1.4 using equation 1.7 



dB 

~dt 



vE w^S^i^^^i B 

c Iq Iq 



(1.12) 



where to and lo are the time and length scales of the plasma. Thus equation 1.4 can be 
reduced to : 

V X B = 47rJ. (1.13) 



Solving equation 1.13 for the current and substituting the result in equation 1.9, we have 



p^- + v.Vv 



-Vp+ —(V xB) xB + pg 



(1.14) 
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Thus, the basic equations of ideal MHD are : 

|^ + V-(pv) = (1.15) 

p( ^ + v-Vvj =-Vp+— (VxB) xB + pg (1.16) 

^=Vx(vxB) (1.17) 

V-B = (1.18) 



Another equation often used by the models is the Lorentz force. The Lorentz force for a 
infinitely conducting plasma can be divided into two terms, magnetic pressure and magnetic 
tension. 

Fl = ^(VxB)xB = -3^B-VB--^VB2 (1.20) 
47r 47r Stt 

where B • VB is the magnetic tension and V-B^ is the magnetic pressure. It is often useful 
to express the Lorentz force in this way because the magnetic pressure is analogous to the 
gas pressure. Another useful plasma parameter is the Alfven speed. 



It is the electrodynamic equivalent of the sound speed in a gas. 
1.3.2 Kinematic Model 

The kinematic model is a simple view of CME propagation and expansion. In the kinematic 
model the CME is considered as defined volume enclosing a fixed mass. Thus, the forces on 
the CME volume can be solved for by integrating over the surface of the volume. Here, we 
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will review two example of kinematic models; melon-seed and melon-seed overpressure. The 
kinematic model offers little insight to the physical properties of the CME since all motion 
is being driven by forces outside of the CME. However, the simplicity of the approach is 
attractive. If this simple approach can describe the expansion and propagation of the CME 
then perhaps the output of the kinematic model could be used in more sophisticated models 
to understand physical parameters of the CME. 

Melon-Seed Model 



In the melon-seed model, Pneuman (1984) describes the CME as a magnetically isolated 
structure in the ambient solar coronal magnetic field. The magnetic structure of the CME 
defines a volume of constant mass. The forces acting on the volume of the CME are the 
ambient magnetic field, the ambient gas pressure and the solar gravity. The Lorentz force 
of the external magnetic field exerts both a magnetic tension and pressure force on the 
CME. The magnetic tension on the CME tethers it to the solar surface. The gradient of 
the ambient gas pressure and the magnetic pressure of the coronal field are pushing the 
CME away. The gradient of the gas pressure is assumed to be insignificant compared to 
the magnetic pressure. Thus the visually descriptive name of the model comes from the 
coronal magnetic field squeezing the CME away from the Sun. 

From the ideal MHD equation, the contributions of the gas pressure and Lorentz force 
can be written as 

F = -Vp + -"-(V X B) X B. (1.22) 

Air 

We can integrate the force expression over the entire volume of the CME structure, thus 
the net force on the center of mass can be expressed in the form of a surface integral 



?2 



i^(Be-cis)Be-(Pe + ^Hs 



(1.23) 



where the subscript e refers to the conditions external to the CME. The external Lorentz 
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Figure 1.8: Geometry of magnetic structure Pneuman (1984). The magnetic field lines 
connected to the photosphere tether the CME to the Sun while the magnetic pressure lifts 
the CME upward. 



force has been expanded using equation 1.20 The magnetic pressure force can be converted 
into a volume integral of the form 



F^= ///-V(g)dy. (1.24) 



If we assume that the ambient field outside the CME can be approximated as constant 
around the CME, then the magnetic pressure force is 



Fm^-V-V{^]. (1.25) 



Thus it is the total volume of the CME, not the specific shape, which is most important in 
determining the net force. 

The derivation of the magnetic tension force is much more involved and does not effect 



the CME expansion. For a magnetic structure like that shown in figure 1.8, the top section 



18 



{R < Ri) of the CME does not contribute to the magnetic tension. The only contribution 
to the magnetic tension is from the field lines less then Ri which are tethering the CME to 
the solar surface. 



Pneuman (1984) goes on to estimate the expansion of the CME with respect to the 
height for a pressure dominated and magnetic dominated interior assuming a force balance 
at the surface of the CME. For the pressure dominated case, the internal gas pressure must 
balance the external magnetic pressure. Thus, we can write the volume of the CME as 



y = Tir^Rcj) = — = i ^ (1.26) 

P P Bi 



where M is the mass of the CME, Vg is the internal sound speed, and (j) is the angular extent 



of the arcade (perpendicular to figure 1.8). If we solve for r and express the magnetic field 



as radial B(, = BqRq/R^ where Bq is the field strength at the solar surface, then the radius 



■-0 

of the CME goes as 



BoRl' 



(1.27) 



Thus for the internally pressure dominated case, the CME expands faster than it propagates 
away from the Sun. For the other limiting case, the magnetically dominated interior, the 
radius goes as 



where Tq is an arbitrary function. From observations of ICMEs, we see that the pressure 



dominated case is an incorrect description of the kinematics. In Lepping et al. (2008), the 
average speed of the MC is 458 km s~^ which is an order of magnitude larger then the 
expansion speed. Thus the approximation that the gas pressure dominates the interior of 



the CME seems unlikely at 1 AU. However, Colaninno and Vourlidas (2006) shows that low 



in the corona some CMEs have bulk speed less then the expansion speed. 
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Melon-Seed Overpressure Expansion Model 



The melon-seed overpressure expansion (MSOE) model (Siscoe et al. 2006), uses a very 



similar derivation to find the propagation and expansion of the CME. In the MSOE model, 



as m 



Pneuman ( 1984 ), the gradient of the ambient magnetic field pushes the CME away from 
the Sun as well as the external gas pressure. The MSOE model describes fast CMEs that 
have been magnetically released from the photosphere and ignores any remaining magnetic 
tethering forces. The model is parametrized by the magnetic force pushing the CME up 
opposed by gravity and drag pulling it down. The radial expansion is parametrized by the 
overpressure interaction of the internal and external media using standard fluid dynamics. 
This model strives to construct a minimal set of equations for which the propagation and 
expansion speeds are the outputs. 

The propagation equation is a force equation with four terms that characterize the forces 
due to magnetic repulsion, inertia, gravity, and drag 

Fnet = Fm-Fg- F^. (1.29) 

The force due to gravity is 



Fg = Mqo-^, (1.30) 



where M is the mass of the CME, R is the height of the CME center, Qq is the gravitational 
acceleration at the solar surface and Rq is the initial height of the CME. The magnetic 
repulsion force is, 

F^ = VolBl^ (1.31) 

where Vol is the volume of the CME and Bo is the magnetic field at the surface. The final 
term of the propagation equation is the drag force due to the solar wind. 



Fd = ^Cd{R)Area{r)pa{v - Vsu,f, (1.32) 
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where Vgw is the solar wind velocity and Cd{R) is the drag coefficient of solar wind. Though 



not denoted in Siscoe et al. (2006), this force also requires a surface area which would be 



the surface area of the CME. 



To calculate the expansion rate of the CME, Siscoe et al. (2006) borrows directly from 



classical hydrodynamics. The magnetic field inside the CME is assumed to be stronger 
then the external field of the ambient corona and causes the CME to expand radially in all 
directions. The CME is magnetically isolated from the surrounding field, thus to conserve 
magnetic flux the magnetic pressure must decrease as 



PB = PBo 



(1.33) 



The expansion is taken directly from classical hydrodynamics, only replacing the adiabatic 
gas pressure with the magnetic overpressure. Thus the CME expansion is. 



r(t) 



dV(t) 



(dr{t)y f To 
\ dt ) ^ \r{t) 



{A 



ho 



Po 



pa 



(1.34) 



where Aho is a hybrid Alfven speed based on the initial overpressure field strength within 
the CME and the initial ambient density, pao outside the CME. 

This formulation assumes that the radius of the CME is smaller than the scale height of 
the density of the ambient corona. In other words, the ambient density can be approximated 
as a constant. Also this formulation assumes that the ambient medium is incompressible. 
While this is not a good approximation, it mostly effects the kinematics of the surrounding 
medium which were are not addressing here. 



With the addition of a solar wind model and a coronal density model, equations 4.13 



and 1.34 form a set of coupled ordinary differential equations. The only input parameters 
that are needed are the external field strength, Ba, internal field strength, Bq, and the 



CME mass. For nominal input parameters, ( Siscoe et al. 2006 ) modeled the radius of the 



CME as it expands out into the interplanetary medium. For three values of the internal 
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Figure 1.9: ICME half-thickness as a function of distance from the Sun calculated for a set 



of initial field strengths ( Siscoe et al. 2006 ) 



magnetic field, the radii obtained (figure 1.9) compare well with measurements of ICME 



half-thicknesses at Earth and Venus ( [Owens et al. 2005, Russell and Mulligan 2002). 

The MSOE model uses the strength of both the internal and external magnetic fields, 



unlike Pneuman (1984) where only the external field is considered. Thus this formulation 
of the kinematic model offers more insight to the physical properties of the CME. However, 
many more inputs are necessary since the coronal density and solar wind drag are included. 



1.3.3 Self-Similar Magnetohydrodynamics Model 

In contrast to the Kinematic models, the self-similar models focus on the magnetic field 
internal to the CME. The solutions derived from self-similar MHD can be applied to a variety 
of magnetic structures. Most offend the force-free solution is applied to the cylindrical of the 
magnetic cloud ( Burlaga et al.|[l981 ). Self-similar MHD models rely on directly solving the 
MHD equation for an azimuthally symmetric magnetic field that expands in a self-similar 
way. Self-similar solutions are often the leading term in an asymptotic expansion of a non- 
self-similar evolution, in a regime far from the initial conditions and far from the infiuence 
of boundary conditions ( Low|1982a ). For this model, we will solve it for two different sets of 
initial conditions. First we will assume a force-free configuration where there are no forces 
acting on the plasma of the CME. This scenario is most applicable to magnetic clouds or 
ICMEs that have propagated away from the Sun into the interplanetary medium. Closer 
to the Sun, these assumptions are less applicable. The more general solution of the model 
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is to assume that there is pressure gradient and gravity acting on the CME. We will derive 
the solution for the self-similar MHD model for both situations. 

The self-similar MHD model assumes a force-free magnetic field, define as 



(V X B) X B = 0. 



(1.35) 



The force-free condition requires that 



(V X B) = aB. 



(1.36) 



where a is some function of position. The self-similar MHD model assumes the case where 
B is symmetric about the azimuthal axis and a is constant. In cylindrical coordinates, a 
force-free magnetic field with no azimuthal dependence can be written as. 



B = z X rBt + V X (z X rBj, 



(1.37) 



the sum of a toroidal and poloidal magnetic field which do not depend on the azimuthal 



angle, Bt{r,z) and Bp{r,z) ( Chandrasekhar 1956). Thus the magnetic field is 



B = -r^Bpf + rBt9 + -^[r^Bp\z 
oz r or 



(1.38) 



and the curl of the magnetic field is 



V X B 



-r—Btf + rA55p0" + - — [r'^Bt]z 
•"Or 



dz 



r or 



(1.39) 



where 



A, 



3 5 d'^ 9^ 
r dr dr"^ dz"^ 



(1.40) 
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For equation 1.36 to be satisfied 



(1.41) 



Thus the magnetic field is force-free if both the toroidal and poloidal fields are scalar func- 
tions. 



Free Expansion : Magnetic Clouds and ICMEs 

To formulate the self-similar MHD solution for a magnetic cloud or ICME, we begin by 
assuming that the influence of solar gravity is negligible. If, as we showed in the previous 
section, the magnetic fleld configuration is force-free and the gas pressure gradient is also 
negligible then we can set all the force terms in the MHD equation of motion (equation 



1.16) to zero, 



/9 ( ^ + V- Vvj = -Vp + — (V X B) X B + pg = 0. 



(1.42) 



If we assume that the velocity of the CME is radial, v = v{r, z,t)r, then equation 1.42 can 
be reduced to 



dv dv Q 
dt dr 



(1.43) 



Equation 1.43 is what Farrugia et al. (1993) calls free expansion. 



If we apply a cylindrically symmetric, force-free magnetic field as defined by equation 



1.38 to the MHD induction equation. 



(9B 

~dt 



V X (v X B), 



(1.44) 



for a radial velocity, the z-component of the induction equation is 



dBt dBt 



dt 



dr. 



(1.45) 
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We now introduce the assumption of self-similarity, which means that the magnetic field is 
a function of only one variable, rather than radius and time, where 



(1.46) 



and ^(t) is a function of time only. We will show that the magnetic field is a function of C 
only more rigorously in the next section. 



Solving equation 1.45 for v, using i?i(C) gives 



r d$ 



(1.47) 



The function of ^{t) can now be solved for by substituting equation 1.47 into the free 



expansion equation [1.43 obtaining 



di \^~dt) ~ \ dt) 



(1.48) 



The general solution to this equation is 



1 9$ 



{t-to)-' 



(1.49) 



where to is a constant of integration. Thus the velocity during free expansion is 



(t + t. 



(1.50) 



If an ICME is moving with a bulk flow speed of U reaches the position of a spacecraft 
at 1 AU at time to, then let t be the time that the spacecraft is inside the ICME. Also let 
To be the radius of the ICME at t = 0. Then as t increases, the distance of the spacecraft 
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from the axis of the ICME is Tq — Ut. The relative speed in the inertial frame is 



U + 



(ro + Ut) 
{t - to) 



(1.51) 



By fitting this equation to the in situ velocity measurements of the magnetic cloud, 



Farrugia et al. (1993) was able to find the propagation time, to, and initial radius, Vo- Using 



these parameters Farrugia et al. (1993) further found that some of the decrease in the 



magnetic field strength seen in the data can be attributed to the expansion of the magnetic 



cloud. Lepping et al. (2008) uses this model to calculate the linear expansion speeds of 



magnetic clouds. 



General Solution 

The general solution of the self-similar MHD model does not require that the gas pressure 
gradient and gravity be negligible. Low published a series of articles outlining a general self- 
similar magnetohydrodynamics (MHD) model ( Low|[l982a|b 1984a|b ). In the first paper of 



the series, Low (1982a) derives an expression for the Lagrangian velocity of the self-similar 



expansion from the MHD equations. We will see that the general solution is similar to 
the free-expansion model. However, more information relating the expansion speed back to 
other physical parameters are found. 

Using spherical coordinates {r,6,(j)), we can express the azimuthally symmetric, force- 
free magnetic field in terms of a toroidal, Bt{r,6,t), and poloidal Bp{r,6,t) fields of the 
form: 

1 



B 



rsmt 



r + Btc, 



r dr 



dr 



(1.52) 
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thus the Lorentz force becomes, 



47r 



(V X B) X B 



or or 



P f QQ r 89 J \r d{r,9) 



(1.53) 



where we have introduced the operator 



sin9 (9/1 d 

C=— + 



dr 



2 89 V sin9 89 



(1.54) 



As before, we assume that the vefocity is only in the radial direction v = v{r,9,t)f. The 

are : 



(r, 9, 0)-components of the MHD equation of motion, 1.16 



/ dv 8v 
'[dt^'^dr 



-1 (^j^ dB, ^ 

47rr'^sin'^9 \ ^ Or 8r 



8p _ GMq 
dr ^ 



(1.55) 



1 (^j^ \dB, ^ j^^ldBA ^ ld_p 
A-Kr'^sirfi9 \ ^ r 89 r 89 j r 89 



(1.56) 



ld{Bt,Bp) 
r d{r, 9) 



(1.57) 



If we apply the magnetic field in equation |1.52| to the MHD induction equation |1.17[ the 
components are : 



(1.58) 



dBt 8 , 



(1.59) 



where the r and ^-components have reduced to the same equation. The MHD conservation 
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of mass equation 3.4 becomes : 



+ — — (rpv) = 



dp 

dt ' dt 



(1.60) 



And finally, assuming that the adiabatic index, 7 = 4/3, the MHD equation of state 1.19 
becomes : 



d f p \ d f p 



dt \ p'^l^ 1 dr\ f?4/3 







(1.61) 



Again we introduce the assumption of self-similarity. We eliminate the r and t depen- 
dence of the equations by defining the control variable 



(1.62) 



where <I>(t) is a function of time only. If we assume that poloidal magnetic field varies with 
time and radial distance through the self-similar control variable C, then : 



Bp{r,e,t)=Bp{C,^ 



(1.63) 



Transforming the variables in equation 1.58, we have : 



r d Id 



0. 



(1.64) 



Solving equation 1.64 for the velocity v, we find : 



r (i$ 



(1.65) 



We have found the same solution as we did for the free expansion, self-similar MHD model 
equation 1.45 However, for this solution we did not assume that the gas pressure gradient 
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and gravity were negligible. This form of the velocity is general to the self-similar model 
assuming a radial velocity and azimuthal symmetry. 



If we make the change of variables in ^-component of the induction equation 1.59 using 
this result, we have : 



r d 1 r d 



(1.66) 



It follows that 



Bt{r,9,t) = -Q{C,e) 



(1.67) 



where Q is any arbitrary function. Equation 1.57 dictates that Bt as a function of r and 6 



can be transformed into a strict function of Bp. Thus we can further restrict Bt as : 



Btir,e,t) = ^QiC,9) = ^Q{Bp[C,e]) 



(1.68) 



Without loss of generality, we can transform any function of r, and t into a function 



of Bp and C (Low 1975). If we transform the pressure, equation 1.56 becomes, 



d^Bp sine d ( 1 dBp\ ^. ,dQ{Bp) 



dC C dd \sine de 



dBy 



(1.69) 



OBr, 



The last term of the equation is the only one explicitly dependent on time. To remove the 



explicit time dependency from equation 1.69 we define : 



p{CBp,t) = —^P{C,Bp) 



(1.70) 
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Making the same variable transform in equation 1.55, we have : 



1 P GMp 



1 



d'^Bj, sinO d 



dC C d9 \sine 89 



dBr, 



(1.71) 



Using equation 1.69 to simphfy this equation we have : 



1 P GMp 



1 



dBr, 



dB„ 



(1.72) 



thus we are left with : 



GMp 

(■2^2 



(1.73) 



To match our other variable, we want p to be dependent on time only through the 
variable Bp. Thus we require : 



a 



_ 



(1.74) 



integrating once 



d^y _r]^-2a 
~dtj ~ ^ 



(1.75) 



The Lagrangian velocity of a particle in the plasma will be : 



dr 



(1.76) 



where Co is a constant. We take Co to the the constant value of C associated with the leading 
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Figure 1.10: Velocity profile for various values of rj and a Low (1982a) 



edge. Thus : 



dr _ /t]^ -2a 



1/2 



(1.77) 



writing the equation with respect to r we have : 



dr ^ ( W ~ 2aCo 

~dt^^°\ T 



1/2 



(1.78) 



Figure 1.10 shows that variety of solutions from setting the integration constants to 
various values. However, for the a = case, the solution simplifies to the free-expansion 
velocity. The general solution of the self-similar MHD model, while more laborious, offers 
a larger array of solutions that are due to the introduction of the gas pressure and gravity. 



1.3.4 Flux Rope Model 



Chen ( 1989 ) explored the toroidal forces as applied to a current loop or flux rope embedded 



in a background plasma such as the solar corona. In this paper, Chen ( 1989 ) only considered 



the MHD aspects of the Lorentz and pressure forces acting on the flux rope. The flux rope 
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consists of a loop wound with a twisted magnetic field. The geometry of the flux rope, figure 
1.11, is defined by the following parameters: Z the height of the center of the loop at the 
apex above the photosphere, R the major radius of the loop, and a the minor radius. The 
flux rope is assumed to extend below the photosphere, though the physics of this region 
are not addressed. The foot-points of the flux rope are separated by the distance 2so and 
are assumed to be fixed in the photosphere. The model is only considered for aspect ratios, 
i?/a, within the range of 5-10. This range of the aspect ratio allows for simplifications of 
the model with respect to the current in the photosphere. 

The geometry of the torus magnetic field gives rise to the specific properties in the 
plasma which do not appear in other magnetic configurations. "Toroidal forces" refers to 
the properties, which may be approximated as a section of a torus. The toroidal forces have 
been fully developed for laboratory plasmas. However, laboratory and solar plasma have 
different boundary conditions. Laboratory plasmas are typically surrounded by a vacuum 
enclosed in a rigid metal wall. Also, magnetic fields are applied externally to balance the 
toroidal force, trapping the plasma. Solar plasma is surrounded by ambient plasma and the 
magnetic field can extend to infinity. 

The flux rope consists of a toroidal Bi and poloidal i?p magnetic field with a poloidal and 
toroidal current density components Jp and Jt-, respectively. The curvature of the magnetic 
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field lines creates a tension force that points inward along the minor radius everywhere. 
In a symmetrically straight cylinder, these forces would cancel out when integrated along 
the cylinder. However, for the torus the top and the bottom field are not symmetric. To 
account for the toroid configuration of the magnetic field, it is necessary to integrate the 
Lorentz force over a toroidal section. The Lorentz force is 

J X B = {JpBt - JtBp)f (1.79) 

where r is the unit vector in the outward minor radial direction. The total force due to 
(—JtBp) is Fp and points in the major radially outward direction, and (JpBt) gives Ft in 
the major radially inward direction. Thus the direction of the Lorentz force depends on the 
difference of the toroidal and poloidal forces. In the solar corona, the pressure force also 
acts on the current loop in the major radial direction. 
The local force acting on a plasma element is given by 



f = Jl(v X B)xB- Vp. (1.80) 



However, to determine the motion of the center of mass of a section of the loop the force 
must be integrated over a given section. To integrate the force over a section of the loop the 
external conditions must be specified: the loop is surrounded by a plasma of pressure pa 
and the magnetic field vanishes at infinity. The integration over the toroidal loop has been 



fully developed by Shafranov (1966). Here we give the equation of motion for the center of 



motion of the top section of the flux rope: 



^ cm 



+ Fg + Fd. (1.81) 



The force, Fr = Md?Z/dt^, is the major radial force per unit length of the loop with a 
mass, M. The magnetic fields, Bt and Bg, are the model ambient fields in the toroidal 
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direction and normal to the plane of the loop, respectively. The quantity (5p is defined by 



/3p 



P-Pa 



(1.82) 



where p is the average internal pressure of the loop, pa is the ambient pressure, and Bp = 
Bp{a) is the poloidal magnetic field at the outer edge of the minor radius (r = a). The 
quantity is the internal induction, characterizing the minor radial current distribution, 
given by 

2 



a?Bl{a) Jo 



r'^Bp{r)dr. 



(1.83) 



The values of vary from for a surface current distribution to 1 /2 for a uniform current 



distribution. It is interesting to note that the form of equation 4.18 was derived assuming 



a circular minor radius of the loop. Also two of the terms of equation 4.18 directly depend 
on the aspect ratio R/a of the toroid. 

In equation 4.18, the terms (Zn(^) — 1 + ^) come from the Jt x Bp force which pushed the 



loop outward. The term {2j^) is from the Jp x Bt force which pushes the loop downward. 

The term (2 ^) is from the external magnetic field which is holding down the loop. 
The direction of the force from the gas pressure /3p can be in either direction. 
The evolution of the minor radius a{t) is given by 



dt"^ Mc^a \ Bj 



l + /3p 



(1.84) 



It is informative to look at the proportionality of equation 4.19 which is 



di2 



OC [B^^-Bl-Bl+STT{p-pa)). 



(1.85) 
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Figure 1.12: Velocity and radius profile Chen (1996). 



Thus, B1 and p cause the minor radius to expand and i?gj, B'^ and pa cause it to contract. 
The effect of the various forces is what we would have expected intuitively. 

If the apex height of a flux rope in equilibrium is perturbed, we can analytically solve 



for the rate of the minor radius expansion in terms of the major radius expansion. Chen 



( 1996 ) gives this proportionality as 



da 


oc 




dZ 


'dt 




dt 



(1.86) 



Thus for the range of aspect ratios considered in the model, the minor radius expansion 
speed is 10 to 20 percent that of the apex expansion speed. 



Figure 1.12 shows the velocity and radius for nominal model values given in Chen ( 1996 ). 
As the flux rope increases in height away from the solar surface, the ambient pressure and 
external toroidal magnetic field will decrease. Applying conservation of flux to the poloidal 
and toroidal magnetic fields, as the minor radius increases, Bt and Bp will decrease. Thus as 
the flux rope moves away from the Sun, the pressure term, 8Tr{p — pa) becomes increasingly 
important. 



1.3.5 Summary of Theoretical Models 

The theoretical models describing the propagation and expansion of CMEs all have ad- 
vantages and disadvantages. The kinetic models are the most intuitive and the easiest to 
understand. However, the melon-seed model seems to overly simplify the problem of CME 



35 



MODELS EXPANSION OBSERVABLES 



Kinetic 

Melon-seed 



MSOE r{t)'^ = -|(^^^)'+ r,To,v,a 



( jro_\ ^^ho~^ao) Po 

\r{t) ) 2 Pa 



Self-Similar MHD 

Free-expansion v = U + v,U,ro,to 

1/2 



General solution ^ = Co (^^^^=^) r, v, (o 

Flux rope 



I da I I I o da 

I dt I \r) \ dt\ dt ' dt 



Table 1.1: Summary of all models reviewed. The observables column lists all parameters 
that can be measured in coronagraphic data. 



propagation and addresses the expansion as an aside. The melon-seed model offers very 
little insight to the physical properties of the CME itself since only the external forces are 
considered. However, as mentioned before, if this simple model can be used to find the 
propagation and expansion of the CME then perhaps these values could be used in a more 
complex model to find other physical parameters of the CME. 

The melon-seed overpressure expansion model uses the simplified geometry of the melon- 
seed model with the addition of several other terms. Clearly, the melon-seed model can be 
improved by some consideration of the physical processes internal to the CME. However, 
it is not clear that the overpressure expansion equation correctly describes these processes. 
The replacement of the adiabatic gas pressure with the magnetic pressure requires more 
explanation and analysis. The addition of the overpressure expansion term which relies 
on the height of CME creates a set of coupled ordinary differential equations which is an 
additional complication to the model since the propagation depends on a large number of 
parameters. 

The self-similar MHD model is far more complicated than the kinetic models. The 
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general solution for the self-similar MHD model incorporates all physical processes of the 
plasma and reduces to the simple solution for the case of free expansion. However, it is 
possible that the added complexity of the self-similar MHD model might make the equations 
too complex to solve. Also the model is still limited by the assumed axially symmetric 
geometry and the force-free magnetic field. 

Last, the flux rope model has many advantages over the other two models. From the 
current STEREO coronagraph images, the geometry that best fits that of the observed 
CMEs is a flux rope. Thus both the geometry and the magnetic fleld configuration of 
the fiux rope model seem to match the CME closely. The fiux rope model has the added 
complication of expansion in both the minor and major radial directions. In the data these 
radial directions are difficult to separate and effect the model in very different ways. Also 
the integration of the force over the flux rope assumes a circular minor radius. As the CME 
moves out into the interplanetary medium the minor radius in no longer circular. How this 
distortion of integration geometry effects the equation of motion has not been studied. 



1.4 Previous Studies 

The expansion speed of CMEs has been measured in coronagraphic images using an image 



processing technique called optical flow. Colaninno and Vourlidas ( 2006 ) used an optical 



flow algorithm to find the expansion speeds of several CMEs observed in SOHO-LASCO. 
The optical fiow algorithm was applied to images taken with SOHO-LASCO C2. The 
authors were restricted to working with data from special observing sequences when the 
image cadence was between 5 to 12 minutes due to the limitations of optical fiow technique. 
Also they were limited to CMEs with height-time speeds less then 300 km s~^. Finally, 
they were only able to measure expansion speed for eight CMEs with large, well-defined 
structure. They found that the sum of the bulk and expansion speeds are all very close 
to the height-time (HT) speed, ranging from within 3% to 18% of the HT speed with an 
average of 9%. These results indicate that the bulk and expansion speeds of CMEs are 
being separated accurately by the optical flow method. 
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Lepping et al. (2008) used in situ magnetic field and plasma data from the WIND 



satellite to calculate the expansion speed of magnetic clouds at 1 AU. Lepping et al. (2008) 



defined the expansion speed as the component of the velocity perpendicular to the MC's 
axes. Using a method derived from the force-free, self-similar expansion of an idealized 
MHD magnetic cloud, the authors calculated the expansion speed of 53 magnetic clouds. 
They found an average speed of 49 km s~^ with a standard deviation of 27 km s^^. The 
estimated error on these speeds is between 15 to 20 km s~^ depending on the MC speed. 



The expansion speed results from Colaninno and Vourlidas (2006) near the Sun have 



an average of 101 km s ^ with a standard deviation of 45 km s ^. These speeds are 



approximately double the expansion speed results from Lepping et al. ( 2008 ) . It seems that 



measurements alone do not complete the picture of CME expansion. Both of these methods 
of measuring CME expansion speeds are limited by geometrical effects. In the coronagraph 
images, all the velocities are projected onto the plan of the sky and are less then the true 
three dimensional velocities. In the derivation of the expansion velocity from the in situ 
data the cross-section of the MC is assumed to be circular even though this geometry is 



unhkely (Lepping et al. 2008). 



1.5 Program of Study 

We begin by selecting CMEs that were well observed in both the remote sensing and in situ 
data between January 2010 and June 2011. In Chapter [2| we discuss the criteria we used for 
selecting CMEs for this study. Of the more than 2,000 entries in the SOHO LASCO CME 
Catalog during this time period, we found only nine CMEs that met our criteria of being 
well observed in remote sensing and in situ data. Also in the Chapter [2| we describe some 
the characteristics of the studied CMEs as observed in the remote sensing and in situ data. 
In Chapter [3| we use the Elliptical and GCS models described in sections 1.2.1 and 1.2.2 
to analyze the remote sensing data. From the GCS model fit, we find the trajectory and 
dimensions of the CME at a range of heights between 3.0 and 211 Rq. The Ellipse model 
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allows us to parameterize any possible elliptical distortion in the cross-section of the CME. 
From the Ellipse model we find the center height and the radius of the CME parallel and 
perpendicular to the CME propagation. At the end of Chapter [3j we compare the results 
of the GCS and Ellipse models. In Chapter [4j we study the propagation and expansion of 
the CME in two ways. First in section 4.1 we analyze the results of Chapter [3] to derived 
the kinematics of the CME. We used the derived kinematics of the CME front to predict 
the arrival time and velocity of the CME at Earth. We compare our predicted values with 
those detected in situ. We also derive the expansion speeds of the CME from the Ellipse 



model results. In section 4.2, we used the kinematics derived in the first half of the chapter 
to calculate the forces on the CME as predicted by the theoretical models introduced in 
section 1.3 In Chapter [5j we use the model described in section 1.2.3 to derive a size and 
orientation of the CMEs from the in situ data. We compare the results from the in situ 
model with the results of the GCS model projected to 1 AU. Finally in Chapter [6j we 
summarize and discuss our results and outline any future work that would lead from this 
study. 
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Chapter 2: Studied CMEs 



In this study, we analyze nine events chosen from all Earth-directed CMEs observed between 
January 2010 and June 2011. The main qualification for the selection of an event was that 
it was well observed in both remote sensing data and in situ data at Earth. In the next 
sections, we will describe the criteria for well observed in each data set. The qualification 
that the CME impact Earth strictly limits the number of possible CMEs. During the 
period of selection, 2,062 CMEs were cataloged in the SOHO LASCO CME Catalog from 
the remote sensing data ( Yashiro et al.|2004 ). During the same time, only 61 possible CMEs 



were identified in the WIND in situ data at Earth (Wu 2011). 

Thus, we began selecting events by attempting to match possible CMEs detected in 
the WIND data to CMEs observed near the Sun in remote sensing data. We did this by 
assuming a three-four day transit time from Sun to Earth and identifying Earth-directed 
CMEs which left the Sun during this time. Matching CMEs seen in remote sensing data 
to CMEs detected in situ is not an easy task even with continuous observation from the 
Sun to 1 AU. Due to projection effects, a CME front may appear to reach a height of 1 
AU several hours before it impacts Earth. Thus if there is more than one CME it can be 
difficult to determine which CME caused the in situ observation. This task is somewhat 
simplified by the relatively low solar activity during the selection time period. Due to the 
low solar activity there are few CMEs during the three-four day window before the in situ 
detections. Also most of the CMEs are isolated in the data. Thus fewer CMEs are obscured 
by another CME and there is less CME-CME interaction. However, the low solar activity 
during the rising phase of the solar cycle also reduces the intensity and size of the events 
during this time limiting the number of CMEs that can be continuously observed in the 
images out to 1 AU. 
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2.1 Selection Criteria for Studied CMEs 



2.1.1 Remote Sensing Criteria 

To select CMEs in the remote sensing data, we first had to identify those that were Earth- 
directed. An Earth-directed CME in the LASCO data wih appear as a partial or full halo. 
In the SOHO LASCO CME Catalog a partial halo is defined as having an apparent angular 
width of more then 120°. With only the LASCO view, determining if a CME is Earth- 
directed is difficult because CMEs that are directed towards and away from Earth have 
the same appearance. Before STEREO, observers had to rely on EUV signatures on the 
solar disk to determine if a CME was directed towards Earth. However, not all CMEs have 



strong or even any EUV signature (Robbrecht et al. 2009) 



With the addition of the STEREO views, it is much easier to determine if the CME 
is Earth-directed. During the period between January 2010 and June 2011, the STEREO 
spacecraft were separated from each other by 132° to 180°. The spacecraft reached opposi- 
tion on 6 February 2011 and afterwards began moving closer to each other on the far side of 
the solar system. On 1 January 2010, the spacecraft were -68°(B) and 64°(A) from Earth. 
On 30 June 2011, they were -93°(B) and 99°(A) from Earth. Given the configuration of the 
STEREO spacecraft during the selection period, an Earth-directed CME would appear as 
a West limb event in STEREO B and an East limb event in STEREO A. Thus, an observer 
can quickly determine if the CME is in the same quadrant of space as the Earth but it is 
still impossible to tell if a CME will produce a measurable signal in the in situ data from the 
remote sensing data. Two of the CMEs we studied were not even identified in the SOHO 



LASCO CME Catalog as partial halos (see table |2.2). 

The first criterion for a CME to be well observed in the remote sensing data is that it 
had to be observed in the 8 telescopes used in this study without a significant period (< 1 
hour) of missing data. Given the nature of space-born observations and the approximately 
three day duration of these CMEs, this first criteria eliminated several events. In all the 
remote sensing data, the CME had to exhibit a bright easily visible structure. Due to the 
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effects of Thomson scattering, the dimmest observations are from the LASCO viewpoint 
for Earth-directed CMEs. Thus the visibihty of the CME in the LASCO coronagraphs is 
usually the limiting case for event selection. The CME had to be easily tracked between 
instruments, thus the structure of the CME had to be visible out to nearly the edge of the 
field of view. Also the CME could not expand outside the upper and lower edges of the 
HI-1 instruments field of view. 



2.1.2 In Situ Criteria 

To identify potential CMEs in the in situ data, we used the automatic detection technique 



of Lepping et al. (2005). The technique was developed to detect potential magnetic clouds 



in the data based on the definition from Burlaga et al. (1981). The technique can also 



identify CMEs in the data. The detection requirements for a magnetic cloud are higher 
then for a CME detection. Thus, we looked at all detections which met the the minimum 
search criteria to find CMEs for this study. The minimum requirements from detection are; 
the proton plasma beta must be low, the field directions must change smoothly and these 
two conditions must persist continuously for a minimum of eight hours. 

Specifically, the CME search was performed on in situ data with a one minute running 
average time step. Periods of low proton plasma beta were identified in the data where 
the average over 25 minutes was < f3p >< 0.3. A smooth change in the direction of the 
magnetic field was evaluated by a quadratic fit of the latitude of the magnetic field, 6b, for 
a running interval of 25 minutes. If the chi-squared of the fit is ^ 500, the magnetic field 
is considered smoothly changing. For a minimum detection, both of these conditions must 



be met continuously for 8 hours. Lepping et al. (2005) has a higher criteria for potential 



magnetic cloud detection. A period of data must meet the minimum criteria as well as 
have a relatively high average magnetic field strength, a low proton thermal velocity and 
the change in the magnetic field latitude must be greater then some lower limit-value. 

We ran the detection software on the WIND data during the search period and found 
61 potential CMEs using the lower search criteria. We searched for CMEs in the remote 
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sensing data that corresponded with these potential in situ detections. We identified nine 
CMEs in the data that match our selection criterion. Of the nine CMEs selected for this 
study, seven also meet the stricter criteria for possible magnetic cloud detections. 



2.2 Description of Studied CMEs 

Throughout the study, we will often refer to the studied CMEs by the date of first appear- 
ance in the LASCO C2 data. We have also assigned a number in chronological order to 
each event. For brevity, we will also refer to the CMEs by these numbers. The date and 



time of the first appearance and event number of each CME are listed in table 2.1 



2.2.1 Remote Sensing Description 

To give a basic description of the studied CMEs, we have chosen some qualitative and 
quantitative characteristics of the CMEs as observed in the remote sensing data. Samples 



of the remote sensing data can be found in Appendix [A} In table 2.1, we have listed the 



widths of each event as seen in the LASCO C2, SECCHI C0R2 A and B images. The 



bottom line of table 2.1 lists the average and standard deviation of each column. All the 



events have a larger width in C2 and similar widths in both C0R2 images as we expect for 



Earth-directed CMEs. In table 2.1 , we also give the speed of the CME front from a linear fit 
of the height-time measurements made for each data sequence. These speeds are projected 
values of the true three-dimensional speed of the CME. We have included the projected 
linear speed for the sake of comparison of these events with previous CME studies when 
only one view point was available. Due to projection effects the same CME can have very 
different speeds from each view point. Several of these events have strong accelerations and 
hence a linear fit is not the best. We will examine the kinematics of the CMEs in depth in 
Chapter |4| 

From the LASCO C2 data, the projected speeds of our CMEs range between 45-669 km 
s^^ with an average and standard deviation of 319 km s^^ ±217. The average speed of our 
CMEs as seen in LASCO C2 is less then the average speed of 437 km s^^ for all halo and 
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Table 2.1: Summary Projected CME Parameters 











Width (deg 







Speed (km s 




CME 


Date 


Time 


C2^ 


C0R2 A 


C0R2 B 


C2^ 


C0R2 A 


C0R2 B 


1 


19-Mar-2010 


10:30 


101 


65 


49 


186 


202 


180 


2 


03-Apr-2010 


10:33 


360 


51 


51 


668 


809 


800 


3 


08-Apr-2010 


01:31 


156 


49 


50 


227 


421 


425 


4 


16-Jun-2010 


14:54 


153 


41 


54 


236 


184 


238 


5 


ll-Sep-2010 


02:00 


236 


70 


78 


469 


392 


406 


6 


26-Oct-2010 


01:36 


83 


71 


73 


214 


189 


260 


7 


15-Feb-2011 


02:24 


360 


90 


91 


669 


751 


773 


8 


25-Mar-2011 


01:25 


85 


50 


48 


45 


100 


142 




25-Mar-2011^ 


14:36 


>191 






119 






9 


Ol-Jun-2011 


18:36 


189 


60 


51 


361 


479 


516 




191±108 


61±15 


61±16 


319±218 392±254 


416±243 



values are taken from the SOHO LASCO CME Catalog 
^ The CME was hsted as two events in the SOHO LASCO CME Catalog. 



partial halo CMEs during any year between 1996 and 2002 in LASCO ( [Yashiro et al.|2004[ ). 
For the same phase of the last solar cycle, 1998, the average speed of all CMEs in LASCO 



was 412 km (Yashiro et al. 2004). This speed is similar to the average speeds as seen 
in the STEREO C0R2 views of our CMEs. Thus the selected CMEs are slow for halo and 
partial halo CMEs but average for all CMEs during this phase of the last solar cycle. 

The average width of the CMEs as seen in STEREO C0R2 is 61° it 16. The qualification 
that the CME not expand out of the HI-1 field of view puts a restriction on the angular 
width of the CME from the SECCHI view point. However, smaller CMEs are usually not 
visible far out into the heliosphere. Thus all the selected CMEs have an angular width 
between 41 and 78 degrees when viewed as limb events. The selected CMEs are larger then 
the average CME as seen in the SOHO LASCO Catalog which detects many small CMEs 
near the Sun ( [Yashiro et al][2004l ). 



In table 2.2 we have placed the events into several descriptive categories. The most 
subjective of theses categories is whether the CME has a fiux rope like structure. A fiux 
rope like structure is defined as fitting the three part structure of an idealized CME. Thus 
the CME must have a rounded leading edge followed by a cavity and bright core. The 
structure of the CME is evaluated in the C0R2 images where it is seen as a limb event. 
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Table 2.2: Categorization of CMEs 



CME 


Flux-rope 


Symmetry 


Halo 


Shock 


Comment 


1 


X 








VpTV ciTciil^ir ^tynctiivp 


2 


X 


X 


Halo 


X 


Bright asymmetric core material 


3 




X 


PH 


X 


Diffuse front 


4 


X 




PH 




Rotating*^ 


5 


X 




PH 


X 


Face-on flux rope structure 


6 


X 


X 






Face-on flux rope structure; diffuse front 


7 




X 


Halo 


X 


Preceding CME; distorted front 


8 


X 


X 


pjja 




Very slow rising; two events in catalog 


9 




X 


PH 


X 


Preceding CME 



Second detection 



^ Vourhdas et al. (2011) and Nieves- Chinchilla et al. (2012) 



Halo events are not seen to have three part structures. From the SOHO LASCO CME 
Catalog, we list whether the CME was observed as a halo or partial halo. Next, we indicate 
whether a white light shock is seen in the data. Shocks can obscure the outline of the CME, 
effecting the model fitting. Shocks are also detected in the in situ data. It is of interest to 
compare shock detection in the two data sets. 

We also evaluated in the C0R2 images the symmetry of the CME between the STEREO 
A and B views. The symmetry gives us an estimate of how close the CME is to the Sun- 
Earth line. The more symmetric a CME is in the STEREO A and B view then presumably 
the closer the CME is to the Sun-Earth line. The symmetry of the CME is also important 
for fitting the structure models to the remote sensing data. If the STEREO A and B views 
are nearly identical then there is little gained from using both view points. However, if the 
asymmetry is large, this might indicate that the CME is distorted. If the CME is distorted, 
it is difficult to fit the idealized shape of the CME models to the data. 

The CMEs chosen through our selection qualifications have a variety of characteristics. 
The projected speeds of the CME are averages for all CMEs. It does not seem that the 
selection qualifications are related to any characteristic other then the size of the events. 
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Table 2.3: Description of CME In Situ 



CME 


Date 


Time 


At 


Rotation 


Shock 


Type 


Dst 


1 


23-Mar-2010 


23:02 


0.67 






1 


1 


2 


05-Apr-2010 


06:43 


1.47 




X 


2 


-76 


3 


ll-Apr-2010 


11:59 


1.30 


By 


X 


2 


-65 


4 


20-Jun-2010 


23:59 


1.70 


Bx 


X 


3 


-6 


5 


14-Sep-2010 


14:24 


1.80 






2 


-22 


6 


31-Oct-2010 


04:48 


1.70 


By 




2 


-7 


7 


18-Feb-2011 


00:00 


3.52 




X 


2 


-25 


8 


29-Mar-2011 


14:38 


3.04 


Bz 


X 


3 


-6 


9 


04-Jun-2011 


19:30 


1.70 


By 


X 


3 


-45 



2.2.2 In Situ Description 



In table 2.3 we have listed some parameters of the studied CMEs as observed in situ. The 
in situ data for each CME are plotted in Appendix |B] For each CME, we have listed its 
arrival time at the WIND spacecraft. This arrival time corresponds to the arrival of the 
jump in density and plasma temperature. The CME arrival time is earlier than the period 
of smoothly varying magnetic field which is used as the arrival time of magnetic clouds. 
Similarly, the end of the CME includes any increase in density following the low density 
cavity of the CME. The duration of the CME in days, At, is given in table |2.3[ We list 
which component of the magnetic field shows a strong rotation and if the CME includes a 
shock. We also list a type for each CME detection. A type 1 CME only meets the minimum 



criteria from Lepping et al. (2005) and a type 2 CME met the higher criteria. A type 3 CME 



is characterized as a MC based on the quality of the fiux rope model from section 1.2.3 
Finally, we list the minimum Dst during the CME passage to indicate the geo-effectiveness 
of each CME. CME 2 with the highest Dst of all the events caused a geo-magnetic storm. 

Of the nine CMEs chosen for this study, three have been studied in previous papers. 
CME 2 on 3-April-2010 generated the first significant geomagnetic storm of solar cycle 24 



and is possibly the most studied CME of the solar cycle to date. Wood et al. (2011) used 



a reconstruction technique to model the CME from Sun to Earth to study the kinematics 



and morphology of the CME. Mostl et al. (2010) studied the in situ signature of the CME. 
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Rollett et al. (2012) and Temmer et al. (2011) studied the kinematics and interaction with 



the solar wind of this CME and two others. Several other papers have been written about 



other aspects of the CME not related to this study (Liu et al. 2011, McComas et al. 2012 



Zuccarello et al. 2012). These works found that the CME is a complex event. Wood et al. 



(2011) found that the orientation of the CME is ambiguous from the remote sensing data. 



Mostl et al. (2010) concluded that the flank of the CME was detected at Earth which kept 



the magnetic cloud from being modeled using standard approaches. CME 4 on 16-June- 



2010 was studied in Vourhdas et al. (2011) and Nieves-Chinchilla et al. (2012). Vourlidas 



et al. (2011 ) found that in the middle corona the CME is rapidly rotating. Nieves-Chinchilla 



et al. (2012) studied the CME in the heliosphere and showed that it continues rotating using 



in situ data from MESSENGER and WIND. In Davis et al. (2011 ), various methods were 
used to predict the arrival time of CME 3 on 8-April-2010 . 

To date, this is the largest study of Earth impacting CMEs using the continuous 



STEREO observation and in situ data. Several studies like Mostl et al. (2010), Nieves- 



Chinchilla et al. (2012), Wood et al. (2011) have studied a single CME from Sun to Earth 



using both remote sensing and in situ data. Rouillard (2011) reviewed these studies which 
incorporated remote sensing and in situ data using a variety of techniques. In this study, we 
apply the same analysis to nine CMEs to derive those parameters which can be compared 
between the remote sensing and in situ data sets. 
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Chapter 3: Analysis of Remote Sensing Data 



To analyze the remote sensing data, we used the Elhptical Model and the GCS Model 



described in sections 1.2.1 and 1.2.2 With these models, we are able to parameterize the 
position and dimension of the CMEs over time from the remote sensing data. In this 
chapter, we present the results of each modeling method. We will further analyze these 
results in Chapter |4] to calculate the kinematics of the CMEs. Also from the remote sensing 
data, we obtained the mass of the CME. We will use the CME mass later to calculate the 
net force on the CME. 

3.1 GCS Model Fit of Remote Sensing Data 

The GCS model is fitted to the remote sensing data by an observer adjusting six parameters 
of the model. The geometry of the GCS model is presented in figure [LB} The parameters of 
the model are longitude, latitude, rotation, ratio, half angle and height. Longitude, latitude 
and rotation are the Euler angles which relate the axes of the model to the Heliocentric 
Earth Ecliptic (HEE) coordinate system. The HEE coordinate system is defined by the x- 
axis along the Sun-Earth line, z-axis is the north pole for the ecliptic of date and the y-axis 
completes the right-handed orthogonal set. The half angle is the angle between the model 
axis and the centers of the legs of the model. The ratio controls the radius of the model. 



The height is the apex height of the front of the model. In table |3.1[ we have listed the 

we 



model parameters for each CME when the model height is ~10 Hq. Also in table 3.1 
have given the width (W) and radius (a) of the model. The geometry of these parameters 
are shown in figure |1.6[ The width and radius are derived from the fit parameters using the 
equations given in |Thernisien| ( |201lD . 



The figures in Appendix [X] show the simultaneously obtained images from the three 
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Table 3.1: GCS Model Parameters at 10 R© 



CME 


Lon 


Lat 


Rot 


Ratio 


Half-Angle 


Width 


Radius 




(deg) 


(deg) 


(deg) 




(deg) 


(Rq) 


(Ro) 


1 


27 


-6 


-35 


0.28 


11 


3.2 


2.2 


2 


6 


-19 


12 


0.37 


25 


4.2 


2.6 


3 


-2 


-3 


-29 


0.19 


33 


4.3 


1.7 


4 


-18 


4 


-33 


0.26 


10 


2.9 


2.0 


5 


-21 


14 


-43 


0.41 


18 


4.0 


2.8 


6 


22 


-28 


-55 


0.27 


30 


4.3 


2.1 


7 


2 


-2 


27 


0.47 


41 


5.1 


3.2 


8 


-27 


3 


-12 


0.31 


41 


4.8 


2.3 


9 


-22 


-5 


-16 


0.39 


30 


4.7 


2.8 



view points as well as the fitted GCS model. In each image, the model is projected onto 
the plane of the image using a grid of points over the surface of the model. For each CME, 



we show the model fit at three heights. The parameters in table 3.1 roughly correspond to 
the first figure for each CME in Appendix [A} 

When the CME is visible in the LASCO data, we use all three view points to make the 
fit. To uniquely fit the GCS model to the remote sensing data, the LASCO view point is 
essential. Since the CMEs are Earth-directed, the geometry of the CMEs as viewed from 



STEREO A and B are similar (table 2.2). The symmetry of the STEREO data set is why 
the LASCO view is of great importance to fit the GCS model. The LASCO view point 
gives us essential information about the orientation and dimensions of the CME which is 
ambiguous from the SECCHI data for Earth-directed CMEs. 

However, given the effects of Thompson scattering and the LASCO field of view, the 
CMEs are not visible in LASCO past ~ 25 R© depending on the density of the front. Thus, 
we must make assumptions about the propagation of the CME to continue fitting the model 
after the CME is no longer visible in the LASCO view. We assume that at large heights the 
CME is expanding self-similarly. This assumption is simple to implement since the model 
expands self-similarly when all parameters except height are held constant. Thus after ~ 
25 Rq, we hold the rotation, ratio and half angle of the model constant. We adjust the 
longitude and latitude to fit any changes in the propagation direction. For most CMEs 
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the model parameters have small variations when fitted using the LASCO view. A notable 
exception is the CME on 16-June-2010 which has a rapid change in the rotation angle in 



the LASCO field of view (Vourlidas et al. 2011). The effects of the rotation on the GCS 



model fit to this CME are discussed in Nieves- Chinchilla et al. (2012). 



We fit the GCS Model at a maximum height of 211 Rq (0.98 AU) for the 3-April-2010 
CME. The average maximum height for all the studied CMEs is 179 (0.83 AU). The 



front height from the GCS model is plotted in figure 4.1 for each event. 

For all the CMEs, the longitude and latitude is within ±30° of the Sun-Earth line. In 
this study we have not associated the CMEs with a solar surface source region. Assuming 
the source regions are similar to the CME trajectory, the locations of the studied CMEs 



are consistent with previous studies of geo-effective CMEs (Zhang et al. 2007). Thus the 



derived three-dimensional trajectory of the CMEs is consistent with all of these CMEs being 
detected in situ at Earth. 



We also note from table|3.1l that the width and radius of the CMEs at 10 R0 are similar. 



The range of widths is 2.9-5.1 R0 and radii is 1.7-3.2 Rq. As we discussed in section 2.2.1 



there was a selection bias for the projected angular width of the CME. This selection bias 
persists for the three-dimensional size of the CMEs as well. 



3.2 Elliptic Model Fit of Remote Sensing Data 



As we discussed in section 1.2, the Ellipse Model allows us to parameterize any possible 
elliptical distortion in the cross-section of the CME. We fit an ellipse to each CME by 
making direct measurements of the CME envelope in the SECCHI remote sensing data. 



Figure 3.1 shows an example of an ellipse fit to SECCHI C0R2 A and B data for the 19- 
March-2010 CME. The direct measurements made by the observer are plotted with plus 
signs. An ellipse fit is plotted over the direct measurements with a solid line. For each fit, 
we get the center of the ellipse, (xc, Uc), the semi-major and semi-minor axis, (a, b) and the 
tilt angle ((/>). The major and minor axes of the ellipse fit are plotted with dashed lines. 
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Figure 3.1: Ellipse fit on the SECCHI C0R2 A (right) and B (left) images froml9-March- 
2010 17:54 UT. The direct measurements (plus signs), ellipse fit (solid line) and major and 
minor axes (dashed lines) have been plotted on the images. 



in figure 3.1 We define the semi- minor axis as the parallel radius and the semi- major axis 
as the perpendicular radius where parallel and perpendicular refer to the direction of the 
CME propagation. 

The results from the ellipse fit at a CME front of 10 R0 from the GCS Model are 



given in table 3.2 in units of Rq. These parameters are from the same images as the GCS 
parameters in table 3.1 To convert the measurements in pixels to distance in units of R©, 
we must know the distance between the measured feature and the observer. Without this 
information, we can only convert the pixel measurements to elongation in degrees. For the 



parameters in table 3.2, we used the CME trajectory from the GCS Model. In table 3.2 



we 



list the conversion factor from Rq to degrees for each spacecraft. The first two parameters 



in table 3.2 are the height of the center of the ellipse fit. 

Even with the correction for the CME trajectory, the ellipse fits in STEREO A and B 
data are sometimes different by as much as 10 R0 for the extreme case of CME 1. There 
are several possible reasons for the differences between the STEREO A and B parameters. 
One possible reason is that the STEREO A and B planes of the sky are different. Thus from 



each view point, we are looking at different projections of the CME. Liewer et al. (2011) 
found a <10° error in the longitude of the leading edge for tie-pointing and triangulation 
technique. They attribute this error to different apparent leading edges. However, for our 
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Table 3.2: Ellipse Model Parameters at 10 R© 



CME 


Center (R©) 


Semi-Major (R©) 


Semi- 


-Minor (R©) 


deg/Ro 




A 


B 


A 


B 


A 


B 


A 


B 


1 


9.3 


8.0 


5.1 


3.6 


2 8 


2 


0.177 


0.263 


2 


9.0 


8.3 


3.5 


3.5 


1.6 


1.8 


0.245 


0.259 


3 


9.7 


10.6 


4.5 


3.8 


1.6 


1.6 


0.261 


0.249 


4 


7.6 


9.8 


2.9 


4.0 


2.4 


2.3 


0.278 


0.205 


5 


7.1 


9.1 


4.4 


6.6 


2.1 


3.1 


0.269 


0.199 


6 


7.9 


7.5 


5.4 


4.8 


2.8 


3.0 


0.243 


0.239 


7 


7.5 


7.7 


4.3 


4.6 


2.1 


2.0 


0.276 


0.258 


8 


9.6 


9.0 


4.1 


3.6 


2.6 


2.3 


0.252 


0.246 


9 


8.0 


8.4 


3.6 


3.0 


2.5 


2.5 


0.250 


0.251 



CMEs, there is no obvious correlation between the differences in the ellipse parameters and 
the spacecraft separation. Also there is no correlation between any of the GCS parameters 
for the orientation of the CME and difference in the ellipse parameters. We will investigate 
the difference in the Ellipse measurements more in the next section by comparing the Ellipse 
and GCS model results. 

We were able to make ellipse fits to the remote sensing data out to a maximum CME 
front heights of 69 R© for the ll-Sep-2010 CME. The average maximum front height for all 
the CMEs is 60 R©. The outer edge of the SECCHI HIl field of view (FOV) is 84 Rq for 
the plane of the sky. Thus we were not able to make ellipse fits over the entire HIl FOV nor 
into the HI2 FOV despite the fronts being visible throughout the HIl and HI2 FOVs. This 
limitation is because the top and bottom edges of the CME diffuse faster than the CME 
front and it becomes difficult to estimate the extent of the CME beyond a certain height. 

3.3 Comparison of the GCS Model Fit to the EHipse 
Model Fit 

We begin by comparing the height of the GCS front with the sum of the parallel radius 
and center height from the Ellipse model. In figure [3^21 the GCS data is plotted with black 
Xs and the Ellipse data for STEREO A and B are plotted with red plus signs and blue 
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squares, respectively. We see that the two values are similar. However, we find that for 
some CMEs (1,4 and 5) the observations from from one view point diverges from the other 
two view pints. In table |2.2[ these CMEs were identified as being asymmetric between 
the STEREO A and B views. Since the direct measurements from one of the STEREO 
views match the GCS fit, we assume that the GCS model was preferentially fitted to this 
view. The preference many not be user error (though that is always a possibility) but 
that one STEREO viewpoint showed a larger deviation from the idealized geometry of the 
model. For example, in figure ??, the CME extends further north in STEREO A than in 
the STEREO B view. The northern extension in the STEREO A data is not fitted by the 
model. Thus the user preferentially fitted the model to the STEREO B data. 

We believe that the discrepancies among the CME measurements is due to intrinsic 
differences in the CME when viewed from different sides. A difference in the CME envelope 
between the STEREO A and B view could be caused by a distortion of the CME by the 
ambient corona. Such a distortion could be caused by a pile up of on one side of the 
CME as it propagates through an non-uniform solar wind. Another possibility is that the 
asymmetry is not a feature of the CME but the projection of another coronal feature such 
streamer which only effects one view point. 



We also compare the radius of the CME derived from the two models. In figure 3.3 



have plotted the mean of the parallel and perpendicular radii and the radius from the GCS 
model using the same convention as figure |3.2[ The radii from the GCS model are similar 
to the mean of the ellipse radii but show the same discrepancies as the front heights. Thus 
the area of cross-section of the two models is also similar. For the theoretical models in 



section 4.2, we will mainly use the GCS values since we have measured the CME with the 
GCS model at greater heights and all the theoretical models assume a circular cross section. 
We use the Ellipse model values when necessary and informative. 



53 




Figure 3.2: Comparison of the GSC front height with the sum of the Ehipse center height 
and parallel radius. From the GCS Model the front heights are plotted with black Xs. From 
the Ellipse Model, we have plotted the center height plus the parallel radius. STEREO A 
and B values are plotted with red plus signs and blue squares, respectively. 
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Figure 3.3: Comparison of GCS radius with the EUipse radii. From the GCS Model the 
radius are plotted in black X. From the Ellipse Model, we have plotted the mean of the 
parallel and perpendicular radii. STEREO A and B values are plotted with red plus signs 
and blue squares, respectively. 



55 



3.4 Measurement of CME Mass 



To calculate the mass of the CME, we used the Thomson scattering geometry from Billings 



(1966) which relates the observed brightness to coronal electron density. This technique 



has been used to measure the mass of CMEs in Poland et al. (1981), Vourlidas et al 



(2000), Vourhdas et al. (2002), Subramanian and Vourlidas (2007), and Colaninno and 



Vourlidas (2009). To convert the number of electrons to mass, we assume that the ejected 



material comprises a mixture of completely ionized hydrogen and 10% helium. To correctly 
calculate the mass of a CME from the Thomson scattering geometry, we must know the 
three dimensional distribution of the CME mass. A first order estimate can be made by 
assuming all the mass is concentrated to the plane of the CME propagation. We used the 
longitude from the CCS Model fit to calculate mass of the CME. 



Figure 3.4 shows the total mass of the studied CMEs versus the height of the GCS model 



front. Colaninno and Vourlidas (2009) found that the CME mass increases with time and 



height until it reaches a constant value above about IORq. To find the final mass of the 
CME, they used the function : 



M{h) = Mc(l - e 



-h/hc 



(3.1) 



where Mc is the final total mass of the event and he is the height where the mass reaches 
63% of its final mass. We have over plotted the fit of this function to the data in figure 



3.4 The results of the fit are given in table 3.3 The final masses are similar to those 



found in Colaninno and Vourlidas (2009). However, the he values (<hc >=1.4) are lower 



then the average <hc >=2.1 found in Colaninno and Vourlidas ( 2009| ). The difference in 
determining deprojected the CME height might account for the difference in he. Here, we 



are using the front height found from the GCS Model fit. Colaninno and Vourlidas (2009) 



used the heights deprojected by the longitude found using their mass technique. 

The results from equation |3. 1| suggest that CMEs are not completely formed in the low 
corona (<1.5 R©). Based on the values of he, the CMEs continue to gain 33% of their mass 
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1 ; 19-Mar-2010 



2 : 3-Apr-2010 
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Figure 3.4: Mass of the studied CMEs as calculated using the Thomson scattering geometry 
and the longitude from the GCS Model fit plotted with squares. The fit of equation to the 
data is plotted with a solid line. 
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Table 3.3: CME Mass 



CME 


he 


log(Mc) 












1 


1.51 


15.3 


2.16 -10^^ 


2 


1.56 


15.6 


4.31 -lO^^ 


3 


1.39 


15.8 


6.13 -10^^ 


4 


1.54 


15.2 


1.77 -IQis 


5 


1.58 


15.8 


5.87 -10^5 


6 


1.42 


15.7 


5.30 -10^5 


7 


1.24 


15.7 


4.92 -10^^ 


8 


1.42 


15.6 


4.22 -10^^ 


9 


1.33 


15.6 


4.18 -lO^^ 



between heights of 1.5 and 10 R©. The increase in mass cannot be completely explained by 
the CME emerging into the SECCHI C0R2 field of view from behind the occulter. This 
increase in mass may be caused by a continued injection of mass from the low corona or 



a pileup of the ambient corona at the CME front. Anzer and Poland (1979) found an in 
flow of mass from lower in the corona for a CME observed in Skylab. Whatever process 
increases the mass, it does not continue past 10 R0. The change in CME mass with height 
will be used when we look at the forces driving the propagation of the CME. 
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Chapter 4: CME Expansion and Propagation 
in the Inner HeUosphere 



4.1 Empirical Results : Elliptical and GCS Model 

In this section we derive the kinematics of the CME based on the measurements made with 
the GCS and Elhptical models. From the GCS model height and time (HT) data, we derive 
the kinematics of the CME front. We use these values to estimate the arrival time and 
velocity of the CME at Earth and compare these results with the in situ data. We also 
use the Ellipse results to find the kinematics of the CME center and the expansion speeds 
parallel and perpendicular to the propagation. 

4.1.1 Propagation 

We used the heights from the GCS Model to determine the three-dimensional kinematics 
of the CME front. The apex of the GCS Model front in Rq is plotted versus elapsed time 



with plus signs in figure 4.1 The elapsed time begins when the CME is visible in LASCO 
C2, C0R2 A and B and the first model fit is made. To analyze the kinematics of the CME 
from the HT data, we must fit an analytical function to these data. We began by fitting 
the data with a first and second order polynomial function. However, we found that a 
single polynomial function could not describe all the HT data. Thus, we have chosen to 
fit the data with multiple polynomial functions. The HT data was fitted using equations 



4.1 - 4.4, The height and velocity of the polynomial equations are continuous by imposing 



the restraints of equations 4.5 - 4.7 This multi-function polynomial fit is used by Wood 



Howard, Plunkett and Socker (2009) and Wood, Howard, Thernisien, Plunkett and Socker 



( 2009 ) to fit the kinematics of two CMEs observed in STEREO from the solar surface to 1 
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AU. 



vot + ho = h 



O.Sait^ + fit + Hi = h 



0.502^ + f2t + H2 = h 



vst + H3 = h 



0<t<ti 

ti <t <t2 
t2 <t <t3 

ts < t < 00 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



where 



Hi = voh + ho 



H2 = 0.5aitl + fit2 + Hi 



H3 = 0.502*1 + /2i3 + H2 



and vi{t = ti) = Vo (4.5) 



and V2{t = t2) = ait2 + fi (4.6) 



and V3 = 02*3 + /2 (4.7) 



The parameters of the multi- function fit are given in table 4.1 For the fit to converge 
some of the parameters must be held constant by the user. The constant parameters for 



each fit are listed in boldface in table 4.1, For all the CMEs the final time, ts, is held 



constant at the time of the last measurement. When fixed by the user, the initial height, 
ho is set at the height of the first measurement. Five of the CMEs do not have an initial 
linear phase at the heights measured. For these CMEs the duration of the initial phase is 
set equal to zero (ti = 0). The total chi-squared of the fit is also given in 4.1 assuming a 
±1 pixel error in the fits for each instrument. 



In figure 4.1 , we have plotted the results for each multi- function fit with the height and 
time data from the GCS Model. Each phase of the fit is plotted in a different color; purple, 
blue, green and red, respectively. The instantaneous velocity is also plotted with a dashed 
line for each CME using the same color scheme. The right axis of the plot gives the scale 
for the velocity data. 
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Figure 4.1: Multi- function fit of time and height data from the GCS Model fit of the studied 
CMEs. The height and time data for each model fit are plotted with plus signs. The multi- 
function fits are plotted in solid lines. Each phase is plotted in a different color; purple, 
blue, green and red, respectively. The velocities of the CME fronts are plotted with dashed 
lines using the same color scheme. 
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Table 4.1: Parameters of Multi- Function Fit 



CME 


ho 


Vq 


ti 


ai 


t2 


a2 


ts 






V ^ } 


(km s~l) 


(hrs) 


(m s"2) 


fhrs) 


(m s"2) 


(hrs) 




1 


3.2 


161.8 


1.9 


8.05 


8.3 


0.35 


77 


2182 


2 


5.0 


916.5 


9.2 


-3.30 


60.0 


0.0 


60 


1921 


3 


3.0 


468.5 


0.0 


2.02 


10.0 


0.0 


71 


23411 


4 


5.7 


193.4 


0.0 


6.47 


7.8 


0.79 


71 


3852 


5 


3.1 


444.3 


0.0 


2.77 


8.1 


-0.68 


58 


7653 


6 


4.3 


214.7 


0.0 


7.20 


10.0 


-0.62 


63 


4780 


7 


4.7 


938.5 


0.0 


-10.77 


11.1 


0.07 


64 


2962 


8 


3.3 


46.7 


8.8 


6.49 


27.2 


-0.16 


107 


14706 


9 


3.0 


1269.5 


0.1 


-19.55 


5.6 


-2.69 


63 


4312 



From the kinematic data, we can immediately see a trend that has been well observed 
in CME propagation ( |Gopalswam^[2006| ). The CME velocities are converging as the CME 
propagates to 1 AU. The range of initial CME velocities is 47-1270 km s~^ with a standard 
deviation of 426 km s~^ while the range of range of final velocities is 321-615 km s~^ with a 
standard deviation of 99 km s~^. The final velocities are listed in table 4.2 Poomvises et al. 



(2010) found similar results. By comparing the velocities in table 4.1 with the projected 



velocities in table 2.1 we see that the initial three-dimensional velocities are larger and have 
a greater range than the projected values. For the LASCO view the initial CME velocities 
range from 45 to 669 km s^^ with a standard deviation of 217 km s^^. Thus the convergence 
of CME velocities from the inner corona to 1 AU is larger by a factor of three when the 
three-dimensional velocities are used. 



Specifically, we see in table 4.1 that the CMEs with the lowest initial velocities (CMEs 
1, 4 and 8) have the largest accelerations in phase 1. While CMEs 2, 7 and 9 have the 
highest initial velocities and are the only ones that decelerate in phase 1. For all the CMEs, 
the acceleration in phase 2 is small, with the exception of CME 9 which has the highest 
initial velocity and decelerated considerably before reaching Earth. CME 9 also has the 
lowest final velocity. As noted in table [2?2] CME 9 is preceded by another CME which many 
account for the unusual kinematics of the event. 

Another extreme case is CME 8. This CME has an initial velocity which is an order 
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of magnitude smaller than the mean of the other CMEs initial velocities. This slow initial 
phase persists for 8 hours before the CME begins to accelerate. By the time the CME 
reaches Earth the finial velocity is similar to the other CMEs. 

Thus all of the CMEs decelerate or accelerate between the Sun and Earth. Only CME 3 
does not have some acceleration in the final phase of the multi- function fits. CME 3 also has 
the largest error for the multi- function fit. The other eight CMEs have a small acceleration 
in the final phase of the fit which persists over an average of 60 hours. An acceleration of 
the CME requires an imbalance of the net force on the CME. We will look at the forces on 
the CME in section [4l 



4.1.2 Predicted CME Values at Earth 

Using the analytical fits to the CME HT data, we can predict the arrival time of the CME 
at Earth. By comparing the predicted arrival time with the in situ data, we can assess 
the accuracy of our three-dimensional trajectories and propagation fits. We can also verify 
the association between the CME observed remotely and in situ in the case of multiple 
events. To estimate the arrival time of the CME, we assume that the CME propagates at 



a constant velocity after the last height measurement, equation |4.4[ Table 4.2 gives the 



predicted arrival times, AT, the difference between the predicted and in situ arrival times 



(table 2.3), and normalized arrival time error, the ratio of the transit time of the CME and 
AT. The average transit time for these CMEs is 3.6 days which is slower then then 2.7 day 
average transit time for CMEs which caused major geomagnetic storms ( Zhang et al.||2003 ) 



but faster then the 4.3 day arrival time for slow CMEs (Gopalswamy et al. 2001). 

For our sample, the maximum discrepancy in the arrival times is -17.6 hours or ±20 % of 
the transit time. This is the largest sample date of arrival times for CMEs using STEREO 



data. Davis et al. (2011) presented the predicted arrival time of CME 3 on 8-April-2010 



from eight techniques. These eight techniques used a range of various lead times and had 



an error of ±12 hours. Our arrival time predictions and those from Davis et al. (2011) are 



not a significant improvement over those made with only the LAS CO data. Gopalswamy 
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Table 4.2: CME Arrival Times and Speeds 



CME 


Arrival 




AT 


AT/Tran 


V3 


L 1 vj vLUi 


Av 




Date 


Time 


(hrs) 


(%) 


(km s^"*^) 


(km s^"*^) 


(km s^^) 


1 


23-Mar-2010 


20:32 


-2.5 


-2 


435 


289 


146 


2 


05-Apr-2010 


17:40 


11.0 


20 


615 


639 


-24 


3 


ll-Apr-2010 


07:46 


-4.2 


-6 


541 


408 


133 


4 


20-Jun-2010 


06:21 


-17.6 


-20 


555 


365 


190 


5 


14-Sep-2010 


22:46 


8.4 


9 


403 


366 


37 


6 


30-Oct-2010 


17:39 


-11.2 


-11 


356 


360 


-4 


7 


18-Feb-2011 


04:32 


4.5 


6 


523 


476 


47 


8 


29-Mar-2011 


14:00 


-0.6 


-1 


430 


342 


88 


9 


05-Jun-2011 


01:04 


5.6 


9 


321 


452 


-131 



et al. (2001) used an empirical model to predict the arrival times of 47 CMEs observed in 



LASCO and found an error range of ±15 hours for 72% of the CMEs. We find that even 
with nearly continuous observations of the CME in the heliosphere and the best available 
models, it is still difficult to predict the arrival time of a CME in situ. 

The only other parameter which can be compared between the remote sensing and in 



situ data without further use of models is the velocity at Earth. In table 4.2, we give 



final predicted velocity, V3 and the velocity in situ at the time of the CME detection. The 
difference between the final and in situ velocity, v^-VinsUu = ^v, ranges from -131 to 190 km 
s~^ with a standard deviation of 99 km s~^. The Av is less then 10 % of the in situ velocity 
for only three CMEs. To see how Av effects the arrival time, in figure |4.2[ we have plotted 
Av versus the normalized arrival time error, AT/Transit. From this plot, we can see that 
the difference in the arrival time and the difference in the final velocity are correlated. If 
the predicted velocity is too fast then the predicted CME arrival time in too early and the 
opposite is also true. 



Based on the results from figure 4.2, we want to test our assumption that the CME 
propagates at a constant velocity after the last height measurement. If the extrapolation 
of the final velocity is inaccurate, we expect the normalized arrival time error to be greater 
with lower maximum heights. In other words, if the extrapolation of a linear final velocity is 
causing the error in the arrival time, we would expect the arrival time error to be small for 
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Figure 4.2: The ratio of the transit time and the error in the arrival time plotted versus the 
difference in the final velocity at Earth. 



CMEs that are measured close to Earth. In figure 4.3, we have plotted the maximum height 
of the GCS Model fit plotted versus normalized arrival time error. The plot shows that error 
in the arrival time is not correlated with the maximum height. Thus our assumption of a 
linear final velocity is not directly effecting the arrival times. 

There are several possible sources of error in the predicted arrival times other than an 
error in the fit of the analytical function to the HT data. One possible source of the error 
is that the predicted arrival times are based on the height of the CME apex. For example, 
if the CME does not impact Earth at the apex the arrival time in situ would be delayed 



assuming a convex front. Mostl and Davies (2012) found that for an assumed circular CME 



front the arrival time delay can be as large as 2 days. However, this effect would only 
produce an early predicted arrival time (AT < 0). Our arrival time errors are both positive 
and negative. Thus this source of error could not apply to all our CMEs. 

Another possible source of error is the fit of the GCS model to the remote sensing data. 
If the longitude of the model fit is incorrect the heights from the fit would also be incorrect. 
In the middle of the HI2 field of view, a 1° change in longitude can cause a change in 
height of ~2 R0 for the same point in the image. The exact error depends on the several 
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Figure 4.3: The maximum height of the GCS Model fit plotted versus the ratio of the transit 
time and the normalized error in the arrival time. 



parameters of the observation. An error in the longitude can make the predicted arrival 
time early or late. Given the distribution of our arrival time error and the sensitively of 
the fit, we assume that this is the most likely source of error. The error in fitting the GCS 
model has not been studied and is beyond the scope of this work. 

4.1.3 Expansion 

In addition to the kinematics of the CME front, we want to analyze the expansion of the 
CME. To isolate the expansion, we decompose the front velocity of the CME into two 



components; the bulk and parallel expansion velocities (Colaninno and Vourlidas 2006). 
The sum of these two velocities is equal to the front velocity. We also want to analyze the 
expansion of the CME perpendicular to the propagation trajectory. We calculate these three 
velocities (bulk, parallel and perpendicular) from the parameters of the Elliptical Model. 
The motion of the center of the ellipse is the bulk velocity. The change in the semi-major 
and semi-minor axes are the perpendicular and parallel velocities, respectively. 

The GCS model has a circular cross-section and thus has only one expansion velocity. 
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The height of the center of the model, R, (figure 



1.6 



point C), is related to the radius, r, by 



r = kR (4.8) 
where k is the user defined GCS model parameter listed in table 3.1 as ratio. If k is constant. 



the model expands in a self-similar way with height. The LASCO field of view is needed for 
K to be well contained. After the CME is no longer visible in the LASCO view, we keep k 
constant. To assess any changes in the bulk and expansion velocities at larger heights, we 
must use the Ellipse Model. 

The results of the ellipse fits are shown in figures \4A 4.6 We have plotted the height of 



the ellipse center versus elapsed time (left column). We have also plotted the parallel and 
perpendicular radii versus elapsed time (right column). The STEREO- A and B values are 
plotted with plus signs and squares, respectively. We have fitted the data with polynomial 
functions. These fits are plotted with solid lines in red and blue for STEREO A and B, 



respectively. The elapsed time along the x-axis is the same as in figure 4.1 

To obtain the bulk velocity, we fit the center height versus time with a second order 
polynomial. Since the ellipse measurements do not cover the range of height that the GCS 



fits do, we found that we need only a single polynomial to fit the data. Table 4.3 gives the 
maximum height of the CME for each ellipse fit, the acceleration and the velocity at the 
maximum height. 

The expansion speeds were calculated by fitting the radii with a first order polynomial. 



The speeds from these fits are give in table 4.4 Also in table 4.4 we list the ratio of the 



parallel and perpendicular speeds. For all the CMEs, the parallel speed is less then half the 
perpendicular velocity with the exception of CME 9 which is about 50 %. These speeds 
confirm what we can see in the images, the CMEs become more elliptical as they propagate 
towards lAU. 



To analyze the elongation of the CME, in figure 4.7, we have plotted the perpendicular 



speeds versus the parallel speeds. STEREO A and B values are plotted with red plus signs 
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Figure 4.4: The height of the elUpse center (right) and the parallel and perpendicular radii 
from the Ellipse model. The STEREO-A and B values are plotted with plus signs and 
squares, respectively. We have fitted the data with quadratic functions. These fits are 
plotted with solid lines in red and blue for STEREO A and B, respectively. All height 
values have been calculated using the longitude found from the GCS Model. 
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Figure 4.5: The height of the elUpse center (right) and the parallel and perpendicular radii 
from the Ellipse model. The STEREO-A and B values are plotted with plus signs and 
squares, respectively. We have fitted the data with quadratic functions. These fits are 
plotted with solid lines in red and blue for STEREO A and B, respectively. All height 
values have been calculated using the longitude found from the GCS Model. 
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Figure 4.6: The height of the elUpse center (right) and the parallel and perpendicular radii 
from the Ellipse model. The STEREO-A and B values are plotted with plus signs and 
squares, respectively. We have fitted the data with quadratic functions. These fits are 
plotted with solid lines in red and blue for STEREO A and B, respectively. All height 
values have been calculated using the longitude found from the GCS Model. 
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Table 4.3: Ellipse Model : Bulk Motion 



CME Max Height Accel. Velocity 





(Ro) 


(m s 




(km 






A 


B 


A 


B 


A 


B 


1 


50.9 


40.5 


2.2 


0.8 


487 


321 


2 


45.9 


50.5 


-13.9 


-3.7 


608 


679 


3 


50.0 


42.8 


-6.3 


-6.8 


339 


379 


4 


42.4 


46.0 


1.7 


1.5 


390 


430 


5 


38.9 


51.3 


-2.1 


-2.6 


340 


401 


6 


40.1 


42.1 


5.0 


2.3 


582 


466 


7 


52.9 


46.8 


2.0 


-5.9 


578 


350 


8 


53.4 


43.0 


3.6 


4.9 


520 


530 


9 


25.3 


26.7 


-30.4 


-34.6 


538 


545 



Table 4.4: Ellipse Model : Expansion Velocities 



CME Perp. Speed Para. Speed Para./Perp. 





(km 


s-i) 


(km 






(%) 




A 


B 


A 


B 


A 


B 


1 


254 


122 


93 


42 


37 


35 


2 


473 


377 


125 


161 


26 


43 


3 


229 


292 


75 


90 


33 


31 


4 


128 


206 


43 


88 


33 


43 


5 


282 


329 


98 


138 


35 


42 


6 


266 


232 


52 


56 


20 


24 


7 


342 


402 


102 


156 


30 


39 


8 


165 


139 


58 


43 


35 


31 


9 


482 


730 


310 


379 


64 


52 
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and blue squares, respectively. We have found that the parallel and perpendicular speeds 
are related by 

Vpara = 0.5 X Vperp - 50. (4.9) 

This result would indicate that the parallel and perpendicular expansions are not indepen- 



dent but are dependent on another parameter of the CME. In figure 4.8, we have plotted the 
parallel and perpendicular speeds versus the bulk acceleration of the CME. Again, STEREO 
A and B values are plotted with red plus signs and blue squares, respectively. The parallel 
and perpendicular speeds are related to the bulk acceleration by 

Vpara = "7.5 X A^ulk " 83 (4.10) 

Vperp = -11.6 X Ai^iu - 250. (4.11) 



Thus for these range of heights where the Ellipse Model was applied, approximately 5 to 50 
Rq, the expansion speeds are constant, the parallel and perpendicular expansion speeds are 
linearly related, and both expansion speeds are linearly related to the the bulk acceleration. 
Since there is no acceleration in the expansion speeds, a result of the ellipse fits is the 
absence of a net force in the interior of the CMEs. The bulk acceleration requires that a 
net force is applied to the exterior of the CME. The result that expansion speeds and the 
bulk acceleration are related, implies that the expansion speeds are also dependent on the 
net external force on the CME. 



4.2 Theoretical Models of CME Propagation and Expansion 

In this section, we use the dimensions and kinematics of the CMEs found from observational 
data with the theoretical models of CME propagation and expansion introduced in section 



1.3 to derive other physical properties of the CMEs. With the theoretical models, we will 



explore the net forces on the CME which drive propagation and expansion. 
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Figure 4.7: From the Ellipse Model, the perpendicular expansion speeds plotted versus the 
parallel expansion speeds. STEREO A values are plotted with red plus signs. STEREO-B 
values are plotted with blue squares. 




Figure 4.8: From the Ellipse Model, the parallel and perpendicular expansion speeds plotted 
versus the bulk acceleration. STEREO A values are plotted with red plus signs. STEREO-B 
values are plotted with blue squares. 
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4.2.1 Kinematic Model 

There are two implementations of the kinematic model, the melon-seed model and melon- 



seed overpressure expansion (MSOE) model, introduced in section 1.3.2 We will use the 
melon-seed model to analyze the expansion of the CME and the MSOE model to analyze 
the external forces on the CME. 

Melon-Seed Model 

The melon-seed model predicts two regimes for the expansion of the CME; the pressure- 



dominated and magnetic-dominated case (Pneuman 1984). The pressured-dominated case 



(equation 1.27) predicts that the CME is expanding faster than the center propagates 



assuming the sound speed is constant (adiabatic expansion). Our measurements (tables 



4.3 and 4.4) show that the pressure dominated case does not apply at these heights in 
the corona. Thus we focus on the magnetic dominated case, which prescribes self-similar 
expansion 

r = r„(|-), (4.12) 

where the radius, r, is dependent on the center height of the CME, R, by some arbitrary 
function Tq- The GCS model expands self-similar ly, where Tq is the user defined parameter 



K, equation 4.8 The parameter k is kept constant after the CME is no longer visible in 
LASCO. 

From the Ellipse Model, we have plotted the ratio of the radii and center height in 



figure |4.9[ The perpendicular data is plotted with squares and the parallel data is plotted 
with plus signs. The STEREO A and B data are plotted in red and blue, respectively. We 
have over plotted the mean of each data set with black lines. The overall trend is that 
CMEs are expanding self-similarly where the ratio of the radii and the center heights are 
constant. Some events, most notably CME 8, seem to show a variation in the parallel 
and perpendicular ratio early in the CME propagation that then becomes constant. We 
investigated this apparent trend and found that variation in the ratio occurs between 10 
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Figure 4.9: The ratio of the parahel and perpendicular radii and center height from the 
Elhpse Model. The perpendicular radius is plotted with squares and the parallel radius is 
plotted with plus signs. The STEREO A and B data are plotted in red a blue, respectively. 
We have over plotted the mean of each data set in black lines. 
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and 15 Rq, the boundary between the C0R2 and HIl instruments. These two instruments 
have different photometeric sensitivities. If the CME has a diffuse front, the outer envelop 
of the CME win appear to increase as it moves into the more sensitive HIl field of view. A 
similar effect has not been observed in limb CMEs as viewed from LASCO C3 where the 
CME is observed in the same instrument until 32 Rq . Thus it is possible that the change in 
the ratios is due to observational effects and is not physical. Other CMEs, such as CME 6 
and CME 9, have variable ratios at later times. However, the variation in the CME 6 data 
is only in the parallel radius and in the CME 9 data is only in the STEREO B data for 
the perpendicular radius. Without a larger sample of CMEs with varying ratios, we cannot 
form any conclusions. Thus we conclude that, most CMEs expand in a self-similar fashion 
in both the parallel and perpendicular directions between 5 and 50 Rq. 



In table 4.5, we have listed the mean of the parallel and perpendicular ratios to center 
height and the k parameter from the CCS Model. The mean and standard derivation of each 
data set is listed in the final line. The perpendicular ratio varies more than the parallel 
ratio. We looked at several parameters of the CME to see if there was any correlation 
with the value of the perpendicular ratio. We did not see any clear trend between the 
perpendicular ratio and the mass, latitude, rotation of the CME or magnetic field of source 
region. The perpendicular ratio may be dependent on structure of the ambient corona. 
Also the K is near the mean of the parallel and perpendicular radii for most of the CMEs. 
This is expected from the comparison of the Ellipse and the CCS fits. 

Melon-Seed Overpressure Expansion Model 



The melon-seed overpressure expansion (MSOE) model (Siscoe et al. 2006) prescribes the 
forces external to the CME. From the MSOE model there are three forces that combine to 
create the net force on the CME and drive its propagation, 



F„.. = M^J = F.„-F,-F.. (4.13) 
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Table 4.5: Ratio of Radius and Center Height 



CME 


Para. 


Perp. 


GCS 


1 


0.27 


0.51 


0.33 


2 


0.21 


0.46 


0.50 


3 


0.19 


0.46 


0.36 


4 


0.24 


0.42 


0.36 


5 


0.30 


0.66 


0.53 


6 


0.29 


0.65 


0.46 


7 


0.27 


0.70 


0.55 


8 


0.29 


0.43 


0.60 


9 


0.33 


0.45 


0.40 



0.27±0.04 0.53±0.11 0.45±0.10 



where is the magnetic force, is the gravitation force, F^ is the drag force, M is the 



mass and ^ is the acceleration of the CME front. 



The gravitational force, Fg equation 1.30, depends on the center height and mass of the 



CME, both measured quantities. The magnetic force, F^ equation 1.31, depends on the 
volume of the CME and the ambient magnetic field. We calculate the volume of the CME 
assuming an ellipsoid shape. 



Vol = -vrr — 
3 2 



(4.14) 



where the radius, r, and width, W, are from the GCS model. To calculate the magnetic 
force, we first need to estimate the ambient magnetic field. We assume that the ambient 
magnetic field is described by the function 



Ba{R) = Bo 



Rq 

R 



(4.15) 



for heights R>5R0 ( Pneuman|1984 ). To find Bq, we solve equation 4.15 where B„(R) is the 
average upstream magnetic field from the in situ data and R is lAU from the in situ data. 



The values of the upstream magnetic field and Bq are listed in table 4.6 The drag force, Fd 



equation 1.32 is dependent on the area of the CME, mass density and drag coefficient of 
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Table 4.6: Ambient Corona 



CME 


Upstream B 


Bo 


Upstream V 




nT 


G 


km s^^ 


1 


3.12 


1.44 


293 


2 


2.40 


1.11 


219 


3 


2.31 


1.07 


235 


4 


3.00 


1.39 


375 


5 


6.63 


3.07 


325 


6 


4.01 


1.86 


369 


7 


3.15 


1.46 


362 


8 


3.69 


1.70 


335 


9 


6.14 


2.84 


365 



the solar wind and difference between the propagation velocity and the solar wind velocity. 



To estimate the density of the solar wind, we use the function from Chen ( 1996 ) 



n{R) 



Rq 

R 



12 



+ 



-^0 

R 



X 10^ + 2.3 X loM % 
R 



(4.16) 



This function gives a density of 5 cm at 1 AU which is consistent with the in situ data. 



To calculate the solar wind velocity, we use the function from Siscoe et al. (2006), 



(4.17) 



where we use the average upstream solar wind velocity from the in situ data for Vswi^ 



AU). The upstream solar wind velocities are listed in table 4.6, We use a drag coefficient 
of 2.5 ( Poomvises||2010 ) . Using these assumed values and estimations of the solar wind, we 
have all the parameters needed to calculate the forces on the CME. 

In figure 4.10[ we plot the logarithm of the absolute value of the forces from the MSOE 
model. We have plotted Fnet, mass times the acceleration of the CME, (solid black), 
(dot-dash black), (solid red) and F^ (dash red) versus front height. The forces plotted 
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in red oppose the propagation of the CME while propels the CME. It is clear from 



figure 4.10 that the drag force completely dominates the other forces by several orders of 
magnitude. Given the values we have used to calculate the forces from the MSOE model, 
the CME would not propagate at these heights. For the forces to equal the net force, the 
ambient magnetic field would need to be 2 to 3 orders of magnitude larger and vary from 



equation 4.15 by an order of magnitude. Similarly, if we were to adjust the drag coefficient 
to balance the forces, it would need to be of the order of 0.001 to 0.01. Thus we must 
presume that there is a fundamental problem in the formulation of the MSOE model. We 
assume that the force imbalance is because there is no consideration in the model of the 
force internal to the CME. The forces internal to the CME could be driving the CME as 
well as the external force. We conclude that considering only the magnetic field external to 
CME in not sufficient to propagate the CME using the observed geometry and kinematics 
from 5 R0 to 1 AU. 

4.2.2 Flux Rope Model 

The flux rope model has the internal magnetic field, gas pressure imbalance and hoop force 
driving the CME as well as the external mag netic field ([Chi^[l989]) . Similar to the MSOE 



model, we can measure or estimate most of the model values. We will use the model to 
calculate the internal magnetic field and compare it to the values detected in situ. Using 
the equation for the force on the center of the CME 



..dv If 



a 



+ Fg + Fd, (4.18) 



the force on of the minor radius 
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Figure 4.10: The forces from the MSOE model F^et mass times the acceleration of the CME 
(solid black), Fm (dot-dash black), F^^ (solid red) and F^ (dash red). 
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equation 1.82 and the definition of the total toroidal current, we find that the internal 
poloidal magnetic field of the CME is 



B, (4.20) 



where 



5 = 2 (^^) Bs (4.22) 
1 AR 

c = 87r{p - Pa) + -aFa -^{Fn-Fg- Fa). (4.23) 



Once we have solved the poloidal magnetic field, we can solve for the internal toroidal 
magnetic field using the equation 

Bt = [aFa - 87r(p - + 5^] . (4.24) 

To use the geometric parameters from the GCS model with the fiux rope model, we 
must reconcile the geometry of the two models. The GCS model creates a croissant-like 
CME where the minor radius varies along the axis. In the flux rope model the minor 
radius is constant along the axis. Also the major radius of the flux rope model is constant 
for all angles. For the GCS model the distance between the center of the model and the 
center of the cavity varies with angle. However, the apex of the two models have the same 
geometry. Thus we use the apex radius and center height of the GCS model for the fiux 



rope parameters, a and Z. We use the equation of Chen (1996) to define the major radius 



R = {Z^ + sl)/2Z (4.25) 
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where Sq is the foot point separation from the GCS model. For all the CMEs the flux rope 
major radius is larger than the width of the GCS model. Thus using the geometry of the 
flux rope model, we are modeling a larger magnetic structure which does not represent the 
size of the density structure seen in the remote sensing data. 



To derive the other parameters in equation 4.20, we begin by solving for the internal 
pressure. Assuming the equation of state is 



where p = nmj is the average internal mass density and 7= 4/3 is the adiabatic index and 
mj is the mass of a proton. Using the geometry of the flux rope to find the density, we can 
also solve for the internal pressure. To find the ambient pressure, we use the equation 

Pa = 2nakTa (4.27) 

where n^ density and Tq temperature of the solar wind. For the ambient temperature, we 
use the equation 

Ta{Z) = To f I] (4.28) 



where To = 2 x 10^ K and a = -0.85 (Chen][l989). We use equation 4.16 for the solar wind 



density. Thus we can solve for the average internal pressure and the ambient pressure. We 
use an internal inductance of = 1.2, which assumes a constant magnetic field distribution. 
We use the same external magnetic field and drag force and gravitational force as for the 
MSOE model. Thus we have all the parameters needed to solve for the internal poloidal 
magnetic field. 



In figure 4.11 we have plotted the magnitude of the internal magnetic field from the 
flux rope model. The dashed line at the bottom is at the value of the maximum magnetic 
field detected at Earth during the CME. For the fiux rope model, there are regions where 

82 



Table 4.7: Flux Rope Internal Magnetic Field 



CME 


H fr 


B fr 


Be 




R,0 


nT 


nT 


1 


118 


14 


7 


2 


210 


100 


14 


3 


195 


82 


12 


4 


142 


10 


8 


5 


100 


22 


11 


6 


44 


133 


12 


7 


189 


7 


18 


8 


94 


24 


15 


9 


166 


10 


29 



the solution for the poloidal magnetic field is imaginary (the square root term in equation 



4.20 is negative). The force balance of the model is very sensitive to the input parameters. 



In table |4.7[ we have listed the minimum internal magnetic field from the flux rope 
model, the height of the minimum field and the maximum magnetic field detected at Earth. 
The results from the flux rope model are very promising. It is possible that by adjusting the 
assumed parameter such as adiabatic index, drag force, and the solar wind parameters, we 
could find values where the flux rope magnetic field would match the measured magnetic 



field at Earth. Kunkel and Chen (2010) were able to reproduce the propagation and mag- 
netic field for a single CME by giving the CME an initial energy profile and adjusting the 



other parameters of the model. Poomvises (2010) was able to reproduce the propagation 
and expansion speeds of four CMEs in the inner heliosphere by adjusting several parameters 



of the model. However, Poomvises (2010) did not have magnetic field measurements for 
their CMEs to compare with the model results. 

Adjusting the parameters of the model may only compensate for the difference in the 
geometry of the flux rope model and the GCS model. The flux rope geometry creates a 
larger CME then we find with the GCS model. Before we can draw any conclusion from the 
flux rope model, we must first investigate how the difference in the observed and flux rope 
geometries effect the outcome of the model. This investigation could be done by deriving 
the flux rope model for the GCS geometry. 
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Figure 4.11: The magnitude of the internal magnetic field from the flux rope model and the 
maximum magnetic field of the CME at earth (dashed). 
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Chapter 5: CME Structure from Remote Sensing and In 

Situ Modeling 



5.1 Modeling of In Situ Data 

Self-Similar Magnetohydrodynamics Model 

We fit the in situ data using the non force free Elliptical Flux Rope (EFR) model described 



in section 1.2.3 The EFR model is similar to the general solution of the self-similar MHD 



model in section 1.3.3 The MHD equations in section 1.3.3 were solved in spherical coordi- 



nates. The EFR model solves the MHD equation in elliptical cylindrical coordinates. The 



relation of elliptical cylindrical coordinates to cartesian is given by equation 1.2 Solving 
the MHD equation in elliptical cylindrical coordinates, allows the model to find solutions 
from the in situ data where the cross-section of the flux rope is elliptical. 

The EFR model is fitted to the in situ data for the regions of the CME where the 
magnetic field is varying smoothly and the proton /3 is low. The boundaries of the EFR fit 
regions are shown in the figures in Appendix |B] with vertical dashed lines. Only three of 
the CMEs (4, 8 and 9) show enough rotation in one of the magnetic field directions to be 



classified as magnetic clouds (table 2.3). We expect the EFR model fit to these CMEs to 
be better than the other CMEs that show little or no rotation in the magnetic field. The in 
situ data used in the model was smoothed to four magnetic field measurements per hour. 

The model parameters that we investigate here are those that describe the size and 
orientation of the CME. The orientation of the CME is defined by the three Euler angles {9, 
(j), ^) that relate the Y, X and Z axis, respectively, of the CME to Geocentric Solar Ecliptic 
(GSE) coordinates. For the remote sensing data, we have been using the Heliocentric Earth 
Ecliptic (HEE) coordinates. The GSE coordinates are centered at Earth and mirror the 

85 



Table 5.1: EFR Model Parameters 



CME 


Lon (p 


Lat e 


Rot i 


Major 


Minor 


Yo 




(deer) 


(dee) 


(deg) 


(AU) 


(AU) 


(AU) 


1 


206 


-44 


108 


0.35 


0.33 


0.34 


2 


119 


-66 


12 


0.30 


0.29 


0.23 


3 


358 


-5 


28 


0.02 


0.01 


0.02 


4 


41 


-14 


174 


0.10 


0.09 


0.05 


5 


19 


-31 


16 


0.36 


0.08 


-0.27 


6 


179 


-69 


145 


0.20 


0.14 


0.03 


7 


63 


-61 


384 


0.37 


0.36 


0.14 


8 


1 


7 


89 


0.04 


0.01 


0.01 


9 


273 


-18 





0.09 


0.08 


-0.08 



HEE coordinates which are centered at the Sun. Also the vector, Yo gives the impact 
parameter, the distance of the spacecraft crossing from the center of the flux rope. The 
semi-minor and semi-major axes define the cross-section of the CME. The results of the 
EFR model fit are listed in table 15. II The best results from the EFR model are for those 
CMEs that have been identified from the in situ data as magnetic clouds, CMEs 4, 8, and 



9. These results are highlighted with boldface text in table 5.1 



5.2 Comparison of the GCS Model with the EFR Model 

Given the different coordinates systems and geometry of the GCS and EFR models it is not 



straight forward to compare the size and orientation results. In figures, 5.1, 5.2 and 5.3 
we have plotted representations of the two models in three planes. The left column is the 
elliptic plane (top view) where the Sun is at (0,0) and Earth is at (1,0). The center column 
is the Sun-Earth line plane (side view) where again the Sun is at (0,0) and Earth is at (1,0). 
The right column is the Earth observer view (head-on view) where the Sun and Earth are 
at (0,0). 

In each view, we have plotted the projection of the elliptical cylinder of the EFR model 
in green. The solid ellipse crosses Earth. The dashed ellipses are it 0.1 AU from the center 
ellipse. The darker ellipse is the same in all three views to give the correct sense of depth. 
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Figure 5.1: Comparisons of the size and orientation of the CMEs 1, 2, ami 3 from the GCS 
Model and the EFR Model. 
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Figure 5.2: Comparisons of the size and orientation of the CMEs 4, 5, and 6 from the GCS 
Model and the EFR Model. 
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Figure 5.3: Comparisons of the size and orientation of the CMEs 7, 8, ami 9 from the GCS 
Model and these EFR Model. 
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Table 5.2: GCS Model Parameters 



CME 


H 


Lon 


Lat 


Rot 


Radius 


Width 




(Rq) 


fdeg) 


fdeff) 


(dee) 


(AU) 


(AU) 


1 


1.66 


26 


-3 


-35 


0.41 


0.60 


2 


1.01 


8 


-2 


12 


0.32 


0.50 


3 


1.10 


-10 


8 


-50 


0.28 


0.49 


4 


1.22 


-17 


6 


-33 


0.32 


0.41 


5 


1.31 


-21 


17 


-43 


0.46 


0.59 


6 


1.33 


13 


-29 


-55 


0.42 


0.62 


7 


1.00 


2 


3 


26 


0.36 


0.51 


8 


1.06 


-20 


3 


-11 


0.40 


0.59 


9 


1.12 


-21 


-3 


-16 


0.29 


0.50 



For the GCS model, the blue outline is the maximum measured height. The black 
outline is the model projected to the height when the front would cross Earth. The dashed 
line is the axis at the center of the cavity. In the left, we have plotted the face-on outline 



of the GCS model, figure 1.6 In the center column, we have plotted the edge-on outline of 



the model. In the right column, we have plotted a top-view profile where the major axis of 



we 



the ellipse is the GCS width and the semi-major axes it the GCS radius. In table 5.2 
list the GCS parameters of the model projected to Earth crossing. 

When comparing the sizes of the CMEs, we expect the results of the EFR model to 
be smaller then the GCS model. From the in situ data in Appendix |B| we can see that 
the density structures have a larger transit time than the magnetic structures. Assuming a 
uniform distribution of the density around the magnetic structure and that both structures 
are traveling at the same speed, we would expect the ratio of the transit times to be the 
same as the ratio of the sizes of the two structures. 

CME 1 : According to the GCS, fit the CME only hits Earth with a glancing blow. 
From the head-on view (right column) the tilt of the CME is different by only 9° comparing 
the two models. Also the radius of the CMEs from both models are different by 0.08 AU. 
However, the impact parameter from the EFR model does not match the GCS model. 
CMEs with large impart parameters are the most difficult to fit with the EFR model since 
the in situ data shows very little rotation in the magnetic field to constrain the fit. 
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CME 2 : The size and impact parameter of the CME is similar from both models. The 
difference in the radius from the two model is 0.02 AU. However, the orientation from the 
two models is very different. From the GCS data, the CME has an East- West orientation 



while the EFR model result has a North-South orientation. Wood et al. (2011) found the 
same difference in orientation for this event between the Sun and Earth. Given these two 
independent results it is possible that the CME has rotated between 25 Rq and lAU. 

CME 3 : There is no similarity between the results of the two models for this CME. For 
this CME, the size of the CME between the two models differs by an order of magnitude. 
The tilt of the two models as seen in the head-on view (right column) is different by 45°. 

CME 4 : The results of the two model are most similar for this CME. The in situ 
detection of this CME is classified as a magnetic cloud. The magnetic structure from the 
EFR model is smaller then the density structure from the GCS model as we would expect. 
Also the orientation of the EFR model in nearly parallel to the front of the CME in the 
ecliptic view (left column). The tilt of the CME is different by 19°. 

CME 5 : From the EFR model the cross section of this CME is highly elliptical. From 
both models, the cross section of the CME is large. However, from the in situ data for 
this CME the density structure is five times larger then the magnetic structure. All other 
parameters are similar between the two models for this CME with the tilt angle which 
differs by 12°. 

CME 6 : This CME is an another glancing blow at Earth. The orientation of the two 
model fits are at very different as well as the impact parameter. The EFR structure is 
orientated North-South while the remote sensing structure is at 55° from the axis. 

CME 7 : This is the only CME where the radius from the EFR model is larger then 
the GCS fit. We would not expect the magnetic structure to be larger then then density 
structure from the in situ data. Also the EFR model intersects Earth at a large angle in 
the first panel while the GCS model has the CME intersecting Earth almost straight on. 

CME 8 : Again for this CME, the EFR radius is an order of magnitude smaller then the 
GCS model while density structure in situ is only 1.4 larger then the magnetic structure. 
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Despite the difference in size, tfie orientation of the CME from both models is similar. The 
difference in the tilt is 3°. This CME was identified in the in situ data as a magnetic cloud. 

CME 9 : This CME was also identified as a magnetic cloud. The orientation of the two 
models is different. However, for this CME, we might be detecting the leg of the CME with 
the EFR model. If we are seeing the leg of the CME then orientation of the EFR model 
pointing back to the Sun is what we would expect. The orientation of the GCS model agrees 
with this interpretation of the two model results. 

Prior to this study, the comparison of the remote sensing density structure to the mag- 
netic structure from the in situ data has not been attempted. Given the different data 
sets and the difference in geometry and coordinate systems of the two models, it is difficult 
to quantify this comparison. The radii from the EFR model vary between 0.01-0.37 AU 
while the GCS radii only vary between 0.28-0.42 AU. With only the results from these two 
models, we cannot know if the magnetic structure of CMEs varies more than the density 
structure or if the results from either model are incorrect. However, looking at the ratio of 
the transit times of the two structures and the ratio of the radii from the two models, we 
suspect that the models are incorrect. 

For the CMEs that arc identified as magnetic clouds (4, 8 and 9 ), the radii from the 
EFR model are less then 0.1 AU. The EFR model is best fitted for these CMEs. Considering 
only these CMEs, we would conclude that the density structures seen in remote sensing data 
are 3-10 times larger then the internal magnetic structure. However, the transit times of 
the spacecraft through the structures does not support this conclusion. The transit times 
of the density structure are only 1.5-5.7 times larger then the magnetic structure. 

The disagreement between these two models has large implications to space weather since 
the orientation of the CME's magnetic field is the other major factor in geo-effectiveness. 
Until the models from these two data sets can be reconciled there in very little progress 
that can be made. Since the in situ data is our only source of direct information of the 
magnetic field and remote sensing data the only means of prediction, it is essential that the 
analysis of these two data sets is reconciled. 
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Chapter 6: Summary and Discussion 



In this study, we analyzed nine CMEs in both remote sensing and in situ data from Sun 
to Earth. By using the remote sensing and in situ data sets, we were able to test many 
of the parameters derived from the models applied to the remote sensing and in situ data. 
From the models of the remote sensing data, we tested predicted values such as arrival 
time, impact velocity and magnetic field at Earth. We also compared the size of the CME 
derived from the in situ data with the density structures seen in remote sensing. By using 
both sets of data for each CME, we can more comprehensively evaluate our analysis of the 
CMEs. We found that there are many avenues for future work based on the results of this 
study. 

The nine CMEs we selected for this study occurred between January 2010 and June 
2011. Of the more than 2,000 CMEs recorded in the SOHO LASCO CME Catalog during 
this time, these were the only nine that met our selection criteria of being Earth impacting 
CMEs that were well observed in the remote sensing data. During the selection period, 
solar activity was increasing from a remarkably low period, the nadir of which was August 
2008. The low level of solar activity is reflected in the number of CMEs we were able to 
track from Sun to Earth in the remote sensing data and the properties of these events. 
From the SOHO LASCO CME Catalog, the projected speeds of our CMEs range between 
45-669 km s^^ with a standard deviation of 217 km s~^. The average speed of these events 



is lower then the average speed of halo and partial halo event as reported in (Yashiro et al. 



2004 ) . We did find a greater range of speeds when we calculated the true three-dimensional 



velocities of the CMEs. The angular widths of the CMEs are between 41 and 78 degrees 
when viewed as limb events. Our selection criteria limited the maximum size of the CMEs 
as viewed from STEREO. By limiting the size of the selected CMEs, we may have limited 



the velocities of the CMEs as well since fast CMEs are often wide (Yashiro et al. 2004). 
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Only three of the CMEs where identified as having magnetic cloud structures detected in 
situ. Overall, the CMEs in this study are unremarkable, average CMEs, with one exception. 
CME 2 on 03-April-2010, caused the first significant geometric storm of the solar cycle 24. 

We began our analysis in Chapter[3]by fitting the remote sensing data with two geometric 
models of the CME structure; the GCS and Ellipse model. The GCS model allowed us to 
estimate the three-dimensional position and size of the CME using the three available views 
of the CME. However, the GCS model assumes a circular cross-section for the CME. From 
the observations, CMEs in the inner heliosphere appear to have elliptical cross-sections. To 
study the cross-section of the CMEs in more detail we used the Ellipse model. The Ellipse 
model is very simple and does not account for projection effects. But this simple approach 
to the data allows us to overcome the limitation of a circular cross section of the GCS model. 
From the Ellipse model, we found the center height and the radius of the CME parallel and 
perpendicular to the propagation. 

In section 3.3, we compared the results of the two models. We found that the sum of the 
center height and parallel radius from the Ellipse model is equal to the front height of the 
GCS model as expected. We also found that the average of the parallel and perpendicular 
radii is similar to the radius of the GCS model. Thus the cross-sectional area of the CMEs 
is similar in both models. This is a helpful result for this study because we found it difficult 
to fit the Ellipse model at heights greater then ~50 Rq. But we were able to fit the the 
GCS model up to a maximum height of 211 R0. To study the kinematics and forces on the 
CME, we wish to use the model parameters at the greatest measured heights. Ultimately, 
our goal is to use the parameters from the remote sensing data in the theoretical models 
of CME propagation and expansion. All of these models assume the CME has a circular 
cross-section. 

We used the height and time measurements from the GCS and Ellipse model in sec- 



tion 4.1, to derive the kinematics of the CME. From the GCS model, we found that the 
kinematics of the CME front could not be well described by a single polynomial equation. 
Thus we used multiple polynomial functions to fit the data. Based on our fit, the CME has 
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different acceleration values during different times in its propagation. All but two of the 
CMEs had an acceleration in the final phase of the its propagation. All the acceleration in 
the final phase are small (< 3 m s^^) but persisted over several hours. We know from the 
in situ data that five of the CMEs do not have the same velocity as the solar wind when 
they cross Earth. Thus these CMEs are not being carried in the solar wind like a leaf in a 
stream. Also the range of CME velocities converge from the Sun to Earth. The range of 
initial CME velocities is 47-1270 km s^^ with a standard deviation of 426 km while the 
range of final velocities at 1 AU is 321-615 km s~^ with a standard deviation of 99 km s~^. 
Thus CMEs must decelerate or accelerate between the Sun and 1 AU ( |Gopalswam^|2006D . 



There does not appear to be a fixed height where the CME reaches a constant velocity. 

An acceleration requires an imbalance of the net forces on the CME. The changes in the 
acceleration of the CME suggests that the net force on the CME is not constant. There 
may be regions where the dominance of individual forces change. It is this force imbalance 
we wish to investigate further with the theoretical models. 



In section 4. 1.2, we used the derived kinematics of the CME front to predict the arrival 
time and velocity of the CME at Earth. We compared our predicted values with those 
detected in situ. The arrival times have a maximum error of ±20% of the total transit time 
and the speeds have an error between -131 to 190 km s~^ with a standard deviation of 
99 km s^^. These arrival times are not better then only using LASCO data (Gopalswamy 



et al. 2001). We found that the error in the arrival time is correlated with the error in the 
final velocity. This result led us to test our use of a constant velocity to project the CME to 
lAU after the final measurement. We did not find that this method was causing the error 
in the arrival times. Delay in the CME arrival time versus the predicted arrival could be 



caused by the shape of the CME (Mostl and Davies 2012). However, our predicted arrival 



times are both earlier and later then the actual arrival time. An error in the arrival time 
can also be caused by a miscalculation of the CME longitude. In the the HI2 field of view 
a 1° change in the longitude changes the height by ~2 Rq for the same position in the 
image. We attribute the overall error in the arrival times to both the estimate of the CME 
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longitude from the GCS model and the fit of the kinematic profile to the height and time 
data. 

To explore the expansion of the CME, we used the results of the Ellipse model. From 



Thernisien (2011), we know that the GCS model is proportional to the center height by 



the factor k which is defined by the user. We fit the parallel and perpendicular radii with 
a linear function. We found that the parallel and perpendicular speeds are related by 
equation |4.9[ Thus the parallel and perpendicular expansions are not independent. This 
result would preclude that the perpendicular expansion is being driven by the solar wind 



while the parallel expansion is being driven by internal forces (Riley and Crooker 2004). 
We also found that the expansion speeds are related to the acceleration of the center height 
of the CME. The acceleration of the center height is related to the net force on the CME. 
Thus both expansion speeds are related to the net force. However, since the expansion is 
at a constant velocity, the net force is not driving the expansion. 



In section 4.2, we used the kinematics derived from the remote sensing data to calculate 



the forces on the CME as predicted by the theoretical models introduced in section 1.3 



The kinematic melon-seed model predicts that the CME will expand self-similarly for the 
magnetically dominated case. We found that most CMEs tend to expand self-similarly. 
There are some deviations from self-similar expansion. However, we cannot determine if 
these deviations are due to observational error or are physical changes in the CMEs. This 
study is focused on Earth impacting CMEs and was not designed to explore if CMEs are 
expanding self-similarly. To properly study self-similar expansion a much larger data set 
would be needed and the CMEs should be measured lower in the corona. From the melon- 
seed model, we can conclude that the interior of the CME is dominated by the magnetic 
field and not gas pressure at heights between 5 and 50 R0. 

From the kinematic MSOE model, we calculated all the forces on the CME as prescribed 
by the model. We estimated the solar wind velocity and magnetic field using the values 
measured in situ upstream of each CME. We found that with the MSOE model, the drag 
force is several orders of magnitude larger then the other forces. The net force on the 
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CME calculated from the measured mass and acceleration cannot be obtained unless the 
magnetic field is 2 to 3 orders of magnitude larger then what we estimated. Similarly, the 
drag coefficient would need to be of the order of 0.001 to 0.01 for the drag force to balance 



the other forces. From comparing our results with those presented in Siscoe et al. (2006), 



they did not use a drag force for most of their examples. Also Siscoe et al. (2006) has a 
maximum radius 0.15 AU for their modeled CMEs at 1 AU where the minimum radius 
we found in the remote sensing data is 0.28 AU. We conclude that there is a fundamental 
problem in the formulation of the MSOE model. We assume that the force imbalance is 
because there is no consideration in the model of the force internal to the CME. 

We then investigated the flux rope model. With the flux rope model, we provide the 
model with all the input values with the exception of the internal magnetic field. We then 
allow the model to solve for this value. From the in situ data we have a measurement of the 
magnetic field to which we can compare the model result. From the fiux rope model we got 



values of the internal magnetic field that are close to those measured in situ (table 4.7). We 
are encouraged by the results of the fiux rope model and see the potential for future work. 

To continue using the fiux rope model with the results of the GCS model we would first 
need to reconcile the geometry of these two models. If we use the dimensions of the GCS 
model apex to generate the fiux rope model geometry, we got a structure that is larger then 
what we observe in the remote sensing data. In the future, we would like to derive the 
equations of the fiux rope model for the geometry of the GCS model. We believe that with 
the correct geometry we would be able to predict the magnitude of the internal magnetic 
field of CMEs from the remote sensing data. 

In Chapter [5j we used the EFR model described in section 1.2.3 to derive a size and 
orientation of the CMEs from the in situ data. We compare the results from the in situ EFR 
model with the results of the GCS model projected to 1 AU. We found that the results of 
the two models do not agree for most events. The EFR model best fits those CMEs that are 
identified as magnetic clouds. For these CMEs, magnetic structures from the EFR models 
are much smaller then the density structure in the remote sensing data. When comparing 
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the sizes of the CMEs we expect the results of the EFR model to be smaller then the GCS 
model. However, we are found that the magnetic structures can be an order of magnitude 
smaller then the density structure. For the CMEs where the sizes from the EFR and GCS 
models are similar, the orientations do not match. The differences between the results of 
these two models is very important. Until the models from these two data sets can be 
reconciled there in very little progress that can be made. 

In conclusion, we have several new results on the kinematics and forces on CMEs as they 
propagate between the Sun and Earth. We found that based on our fit of the height and time 
data, CMEs have different acceleration values during different times in their propagation. 
Also, there is not a fixed height where the CME reaches a constant velocity. From the 
kinematic results, we predicted the arrival time of the CME at Earth and compared this 
arrival time to the detection of the CME in situ. The predicted arrival times have a 
maximum error of ±20% of the total transit time and the speeds have an error between 
-131 to 190 km s~^ with a standard deviation of 99 km s~^. These arrival times are not 



better then only using LASCO data (Gopalswamy et al. 2001). Even with remote sensing 
observations of the inner heliosphere there is still a large uncertainty in the trajectory of the 
CME. We also compared the structure of the CME as predicted from the remote sensing 
and in situ data at Earth. We compare the size and orientation results from the in situ 
EFR model with the results of the GCS model projected to 1 AU. We found that the 
results of the two models do not agree for most events. We also found the kinematics of 
the CME expansion parallel and perpendicular to the propagation of the CME. We found 
that the parallel and perpendicular expansions are not independent. We also found that 
the expansion speeds are related to the acceleration of the center height of the CME. 

We then used the derived kinematics of the CME with three theoretical models of CME 
propagation and expansion. From the melon-seed model, we found that most CMEs tend 
to expand self-similar ly. Also, we can conclude that the interior of the CME is dominated 
by the magnetic field and not gas pressure at heights between 5 and 50 R©. With the 
MSOE model, we found that the drag force is several orders of magnitude larger then the 
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other external forces. The net force on the CME calculated from the measured mass and 
acceleration cannot be obtained with this model unless the magnetic field is 2 to 3 orders 
of magnitude larger then what we estimated or the drag coefficient would need to be of 
the order of 0.001 to 0.01. Our most promising results are from using the kinematics and 
dimensions of the CME with the flux rope model. From the flux rope model, we calculated 



values of the internal magnetic field that are close to those measured in situ (table 4.7). We 
believe that with the correct geometry we would be able to predict the magnitude of the 
internal magnetic field of CMEs from the remote sensing data. 
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Appendix A: Remote Sensing Data 



In this Appendix, we give a sample of the remote sensing data used in the study. For 
each of the nine CMEs in the of the study, we have included the remote sensing data for 
three heights during the CME propagation. In the caption for each figure we give the CME 
number, the time of the remote sensing observation and the height of the GCS model fit. All 
of the figures in the appendix have the same layout. The left panel is data from STEREO 
A. The center panel is data from LASCO. The right panel is data from STEREO B. The 
data in the top and bottom panel are the same. The images in the bottom panel have been 
over plotted with the GCS model. The GCS model is represented by a grid of points on 
the surface of the model. The first figure for each CME has data from STEREO C0R2 and 
LASCO C2. The next two figures have data from STEREO HIl and LASCO C3. 

NOTE : This appendix has been reduced to one figure per event to reduce 
the file size. 
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Figure A.l: CME 1 19-Mar-2010 18:39 UT height = 10 




Figure A. 2: CME 2 03-Apr-2010 10:39 UT height = 6.5 




Figure A.3: CME 3 08-Apr-2010 06:54 UT height = 10 



Figure A. 4: CME 4 16-Jun-2010 19:24 UT height 




Figure A. 5 



CME 5 ll-Sep-2010 05:24 UT height = 10 



Figure A.6: CMR 6 26-Oct-2010 11:54 UT height = 10 
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Figure A.7: CME 7 15-Feb-2011 02:54 UT height = 10 



Figure A. 8: CME 8 25-Mar-2011 11:54 UT height = 10 



Figure A. 9: CME 9 Ol-Jun-2011 20:39 UT height = 10 



Appendix B: In Situ Data 



In this Appendix, we give the in situ data used in this study. In the caption for each figure, 
we give the CME number and date of the in situ detection. All the data is from the WIND 

spacecraft and has been smooth to four points per hour. The hght blue block demarcates 
the boundaries of the CME. The vertical dashed lines indicate the region of the data that 
was used for the EFR fit. 
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Figure B.2: CME 2 05-April-2010 
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Figure B.3: CME 3 ll-April-2010 
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Figure B.4: CME 4 20-June-2010 



114 



^ 1.0!<10' 




^ 1.0 

^ OB 

c 0-6 

£ 0,4 

2 0.2 

°- 0.0 

.-^ SOO 
450 
f 400 

300 



> ^ 150 



< 

E 100 
" ■5' 50 



257.5 



258.0 



25A.5 
DOV Ot £010 



25S,0 



£59.5 



Figure B.5: CME 5 15-September-2010 
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Figure B.7: CME 7 18-Feburary-2011 
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Figure B.9: CME 9 04-June-2011 
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